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PREFACE 


This  study  covers  the  analysis  and  experimental 
investigation  of  the  effect  of  a  hydraulically- supported, 
pivoted-shoe,  journal  bearing  on  the  attenuation  of  noise 
originating  from  rotor  unbalance. 

This  work  is  in  fulfillment  of  Bureau  of  Ships 
Contract  No.  NObs-86914  and  it  has  been  assigned  Project 
Serial  No.  SF  013-11-05,  Task  3679. 

The  final  report  includes  three  volumes.  Volume  3 
being  herewith  included. 


Volume  1  -  Spring  and  Damping  Coefficients  for 

the  Tilting-Pad  Jounral  Bearing. 

Volume  2 

Part  I:  Attenuation  of  Rotor  Unbalance  Forces 

by  Flexible  Bearing  Supports 

Part  II:  Unbalance  Response  of  a  Uniform  Elastic 
Rotor  Supported  in  Damped  Flexible 
Bearings 


Volume  3 

Part  I:  A  General  Computer  Program  for  Unbalance 

Response  of  a  Rotor  in  Fluid  Film 
Journal  Bearings 

Part  II:  Experimental  Investigation  of  Hydraulic- 
Supports 


Mechanical  Technology  Incorporated  was  primarily 
responsible  for  the  analytical  portion  of  this  study,  while 
Westinghouse  Electric  Corporation  designed  and  conducted 
the  experimental  test. 
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Part  I 


A  General  Computer  Program 
for  Unbalance  Response  of  a  Rotor 
in  Fluid  Film  Jounral  Bearings 


by 


Jorgen  Lund 


Mechanical  Technology  Incorporated 
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1NTR0DUCTT.0N 

A  rotor  supported  in  fluid  film  journal  bearings  is  a  complex  dynamical  system 
and  exhibits  a  variety  of  characteristics:  critical  speeds,  instabilities,  un¬ 
balance  vibrations,  etc.  In  many  applications  it  is  a  critical  member  ard  any 
design  procedure  requires,  as  a  minimum,  the  calculation  of  the  critical  speeds 
of  the  rotor.  Other  investigations  may  include  a  check  on  the  stability  of  the 
rotor  (oil  whip,  fractional  frequency  whirl) or  calculating  the  rotor  amplitude 
caused  by  an  external  excitation  (e.g.  shock  loading).  In  the  present  case  the 
concern  is  with  the  response  of  the  rotoi  to  unbalance  forces  and  to  determine 
both  the  whirl  amplitudes  and  the  forces  transmitted  to  the  foundation. 

The  fluid  film  journal  bearings  play  a  very  important  role  in  the  dynamics  of 
the  rotor.  They  are  normally  the  predominant  source  of  damping  such  that  with¬ 
out  this  source  it  would  be  impossible  to  run  the  rotor  through  any  of  its 
critical  speeds.  Secondly,  the  bearing  film  is  flexible  and  thereby  it  may  lower 
the  critical  speeds  drastically  (e.g.  for  the  first  critical  speed  the  reduction 

can  be  30  to  40  per  cent,  or  even  more).  The  fluid  film  flexibility  also  causes 

the  bearing  to  act  as  a  vibration  isolator,  attenuating  the  dynamical  forces 
transmitted  to  the  pedestals.  Hence,  in  any  comprehensive  rotor  response  cal¬ 
culation  it  is  necessary  to  have  a  method  available  which  treats  the  dynamical 
bearing  forces  as  accurately  as  possible.  It  is  the  purpose  of  this  report  to 
present  such  a  calculation  method  and  to  describe  a  computer  program  for  carry¬ 
ing  out  the  numerical  computations. 

The  computer  program  is  very  general.  It  calculates  the  rotor  whirl  amplitude 

and  the  force  transmitted  to  the  base  due  to  a  given  rotor  unbalance.  The  rotor 

is  flexible  and  may  have  any  arbitrary  geometry.  Also,  there  can  be  splined 
couplings  in  the  rotor  and  several  bearings.  The  bearing  pedestals  can  be 
assigned  both  flexibility  and  damping.  Since  the  bearing  film  forces  are  not 
the  same  in  all  directions  the  whirl  motion  of  the  rotor  is  treated  as  two- 
dimensional  such  that  it  becomes  an  orbit  around  the  equilibrium  position.  The 
orbit  is  elliptical  and  its  dimensions  and  orientation  vary  along  the  length  of 
the  rotor.  The  computer  program  calculates  the  whirl  orbits  for  a  number  of 
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points  along  the  rotor  and  gives  also  the  components  of  the  force  transmitted 
to  the  foundations. 

The  report  sets  forth  the  analysis  for  performing  the  calculations  and  describes 
the  computer  program  based  on  the  analysis.  Detailed  instructions  for  preparing 
the  computer  input  and  interpretating  the  output  are  given. 
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DISCUSSION 
a .  General 

The  rotor  analysis  derives  in  its  principle  from  the  Myklestad-Prohl  method 
(Ref.  1,2,3).  However,  in  its  original  form  the  Myklestad-Prohl 
method  is  set  up  only  for  calculating  the  critical  speeds  of 
flexible  rotor  in  flexible  bearings  and  it  treats  the  rotor  motion  as  a 
transverse  vibration  of  a  beam.  In  the  present  case  the  motion  is  treated  as 
two-dimensional,  damping  is  included  in  the  bearings  in  addition  to  stiffness 
and  the  rotor  response  is  calculated  at  any  speed,  not  just  the  mode  shape 
at  the  critical  speeds.  Furthermore,  the  effect  of  gyroscopic  moments  is 
included. 

In  general  a  rotors  cross -sectional  dimensions  and  its  mass  distribution 
varies  along  the  length  of  the  rotor.  Thus,  for  calculation  purposes  it  is 
convenient  to  break  the  rotor  up  into  short  sections,  each  section  having 
a  constant  cross-section.  Furthermore,  when  there  are  many  sections  the  mass 
of  each  section  can  be  divided  into  two  parts  and  lumped  at  the  end  points 
of  the  section.  Concentrated  masses  like  wheels,  impellers,  etc.  can  be  made 
to  coincide  with  an  end  point  of  a  section.  In  this  way  the  rotor  is  replaced 
by  an  idealized  model  consisting  of  a  number  of  mass  points  connected  by 
weightless,  flexible  bars.  The  model  can  be  brought  as  close  to  the  actual 
rotor  as  desired  by  making  the  subdivisions  small  but  in  practice  only  a 
limited  number  of  divisions  is  needed  to  obtain  a  very  good  accuracy. 

Since  the  bearing  film  properties  to  a  large  extent  control  the  whirl  motion 
of  the  rotor,  it  is  necessary  to  represent  the  dynamical  bearing  film  forces 
as  accurately  as  possible.  The  method  of  representation  is  based  on  the 
assumption  that  the  whirl  amplitude  is  small  compared  to  the  bearing  clearance 
such  that  the  dynamical  forces  can  be  replaced  by  their  gradients  around  the 
steady  state  journal  center  position.  In  this  way  the  dynamical  forces  be¬ 
come  proportional  to  the  whirl  amplitude  and  the  corresponding  velocity,  and 
the  factors  of  proportionality  are  called  spring  and  damping  coefficients. 

They  differ  from  conventional  mechanical  spring-dashpot  systems  by  also 
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c ontaining  cross-coupling  terms  in  addition  to  direct-coupling  terms,  i.e.  the 
dynamical  force  in  a  given  direction  (say  the  x-direction)  is  not  only  pro¬ 
portional  to  the  amplitude  and  velocity  components  in  that  direction  but  is  also 
proportional  to  the  amplitude  and  velocity  components  in  the  mutually  perpen¬ 
dicular  direction  (i.e. the  y-direction) .  Hence,  in  an  arbitrary  reference 
coordinate  system  with  x  and  y-axis  the  two  dynamical  force  components  can 
be  expressed  by: 


Fx  -  -  ij  Cxy  ij 

^  *  ~ Kfx  x  ~  ~  kijij  tj  ~  Cijif  ^ 


where  x  and  y  are  the  amplitude  components,  x  and  y  are  the  velocity  components, 

K  and  K  are  the  direct  coupling  spring  coefficients,  C  and  C  are  the 
xx  yy  xx  yy 

direct-coupling  damping  coefficients,  K  and  K  are  the  cross -coupling  spring 

xy  yx 

coefficients,  and  C  and  C  are  the  cross -coupling  damping  coefficients.  These 
xy  yx 

8  coefficients  are  functions  of  the  bearing  Sommerfeld  number  defined  through 
the  rotor  speed,  the  steady  state  bearing  reaction,  the  lubricant,  viscosity  and 
the  bearing  dimensions  (for  gas  bearings  the  coefficients  are  functions  of  the 
compressibility  number  and  the  bearing  eccentricity  ratio).  Thus,  the  coeffi¬ 
cients  va.y  with  speed.  A  method  for  calculating  the  coefficients  is  given  in 
Refs.  4  and  9  and  values  of  the  coefficients  for  several  bearing  types  may  be 
found  in  Refs.  4, 5, 6, 7  and  8. 


Frequently  the  pedestals,  on  which  the  bearings  are  mounted,  are  as  flexible  as 
the  bearing  film.  In  such  cases  the  pedestal  stiffness  must  be  included  in  the 
calculations.  For  completeness  the  analysis  allows  for  both  stiffness,  damping 
and  inertia  in  the  pedestals.  Furthermore,  as  the  rotor  bends  under  the  in¬ 
fluence  of  the  unbalance  forces  the  journals  become  cocked  in  their  bearings. 

The  fluid  film  resist  the  tilting  and  this  can  be  expressed  by  a  set  of  8  spring 
and  damping  coefficients  in  analogy  to  the  previously  discussed  coefficients.  The 
analysis  includes  this  effect,  both  in  the  bearings  and  in  the  pedestals.  The 
resistance  to  tilt  normally  affects  the  rotor  motion  only  at  speeds  above  the 
second  or  third  critical  speed  but  if  the  pedestals  are  made  soft  for  alignment 
purposes  resonance  conditions  may  exist  which  can  only  be  explored  if  the  effect 
of  tilt  is  included. 
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Occasionally  the  rotor  is  not  a  single  member  but  consist  of  several  rotors 
connected  by  splined  couplings  (e.g.  a  turbine-generator  set  connected  by 
a  splined  coupling).  The  analysis  allows  for  including  splined  couplings 
anywhere  in  the  rotor  and  assumes  that  no  bending  moment  is  transfered 
through  the  coupling. 

The  whirling  motion  of  the  rotor  is  generated  by  unbalances  built  into  the 
rotor.  In  general  the  unbalance  varies  in  magnitude  and  circumferential 
location  along  the  rotor  such  that  under  speed  the  unbala  ce  forces  may 
bend  the  rotor  into  complicated  shapes  (e.g.  resembling  a  "cork-screw"). 

The  bend  rotor  whirls  around  its  steady  state  position  (i.e.  the  position 
the  rotor  would  occupy  if  there  were  no  unbalance  forces)with  each  point 
of  the  rotor  axis  describing  an  elliptical  path.  The  dimensions  and  orient¬ 
ation  of  the  ellipse  varies  along  the  length  of  the  rotor. 

If  the  rotor  runs  at  high  speed  and  has  large  discs  (e.g.  turbine  wheels,  etc) 
mounted  on  the  shaft  the  gyroscopic  moment  becomes  important,  especially  if  a 
wheel  is  overhung  at  one  end  of  the  rotor.  The  gyroscopic  moment  is  pro¬ 
portional  to  the  mass-moment  of  inertia  of  the  wheel,  the  square  of  the  speed 
and  the  deflection  angle  of  the  rotor.  If  the  rotor  motion  is  considered  as 
a  transverse  vibration  of  a  beam  (i.e.  the  whirl  orbit  is  a  straight  line) 
the  gyroscopic  moment  tends  to  "soften"  the  rotor  and  lower  the  critical  speed. 
On  the  other  hand,  if  the  bearing  spring  and  damping  coefficients  are  the  same 
in  the  vertical  and  the  horizontal  direction  the  rotor  whirl  orbit  becomes  a 
circle  and  the  gyroscopic  moment  stiffens  the  rotor.  Actually,  the  whirl 
orbit  is  elliptical,  i.e.  somewhere  between  a  straight  line  and  a  circle,  and 
the  effect,  of  the  gyroscopic  moment  can  only  be  assessed  by  performing  the 
complete  rotor  analysis.  It  is  a  non-linear  effect  since  it  depends  on  the 
dimensions  of  the  elliptical  whirl  orbit.  In  the  present  analysis  the  gyro¬ 
scopic  moment  is  taken  into  account  and  is  calculated  by  an  iteration  procedure. 

b .  Special  Considerations  in  Performing  the  Numerical  Calculations 

The  greatest  difficulty  encountered  in  performing  the  numerical  calculations 
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is  the  magnitude  of  the  numbers  and  the  loss  of  significant  figures.  These 
difficulties  become  pronounced  when:  a)  there  is  an  excessive  number  of 
rotor  mass  stations,  b)  the  rotor  is  very  stiff  and,  c)  the  bearings  are 
very  stiff.  There  is  no  universal  remedy  for  the  problem  but  if  trouble 
arises  two  possibilities  may  be  tried:  a)  reduce  the  number  of  rotor 
stations  to  the  essential  minimum  and,  b)  apply  a  scale  factor. 

Let  the  scale  factor  be  d.  .  Then: 
multiply  the  speed  by  Oi 

multiply  (El)  by  ot  (e.g.  multiply  E  by  U  ) 

multiply  the  bearing  spring  and  damping  coefficients  by 

2 

(i.e.  multiply  K  ,<oC  ,M  ,ci>D  etc.  by  o(  ) 

XX  XX  XX  XX 

multiply  the  pedestal  stiffness  by  d*  and  the  pedestal 
damping  coefficients  by 

(i.e.  multiply  Jf*  And  &  by  d  f  ^and  dy  by  U  ) 

leave  the  rotor  masses,  the  rotor  length,  the  pedestal 
masses  and  the  unbalance  unchanged. 

The  the  numerical  results  will  give  the  amplitude  unchanged  whereas  the 
bending  moment  and  the  transmitted  force  must  be  divided  by  d?  to  obtain 
the  actual  values. 


c .  Analysis  and  Dimensionless  Equations 

Referring  to  the  sign  convention  given  in  Fig.  2  and  considering  first  a  con¬ 
tinuous  rotor  the  three  basic  equations  for  determining  the  rotor  motion  are, 

(1-a)  Force  balance  for  a  shaft  increment,  dz  :  s^Aco(x+e) 

(2-a)  Moment  balance  for  a  shaft  increment , dz  ;  f?  *  v+«VV)£ 

die 

(3-a)  Shaft  deflection  :  M  =  dz* 


X  -  amplitude  in  vertical  direction,  inch 

ij  -  amplitude  in  horizontal  direction,  inch 

Z  -  coordinate  along  the  rotor  length,  inch 

t  -  eccentricity  between  mass  center  and  shaft  center,  inch 

2 

A*  cross-sectional  area  of  shaft,  in 


-7- 


I- 

F- 

f 

(ip-lr)  - 
gj  - 

M- 

V- 


cross-section  moment  of  inertia,  in 

2 

Youngs  modulus,  lbs/in 
mass  density,  lbs.sec^/in^ 

mass  moment  of  inertia  per  unit  length,  which  is 

2 

effective  in  gyroscopic  moment,  lbs. sec 
angular  speed,  radians/sec 
bending  moment,  lbs. in 
shear  force,  lbs. 


These  three  equations  may  be  combined  to  give  the  familiar  4th-order 
differential  equation  governing  the  unbalance  vibrations  of  a  rotor: 

(A-a)  -fAtoHx*') 

(see  Ref.  3,  page  330) 

For  a  circular  whirl  orbit: 

(iP-iT)=$r 

For  a  straight  line  orbit: 

For  an  elliptical  whirl  orbit,  see  Eq.(28)  and  (29)  in  this  report. 


At  the  bearings  there  is  an  abrupt  change  in  the  shear  force  and  the  bending 
moment  due  to  the  bearing  reactions.  Let  the  bearing  be  at  Z  =  Z©  .  Then: 

(5-a)  V2wZt+ 

(6 -a)  MZb2+  -  -  ( Mx*  +  ico  Dxx)^- ( Mx> ■+  (Co  Dy) 

where  K  ,  C  ,  M  ,D  etc.  are  the  bearing  spring  and  damping  coefficients. 

XX  XX  XX  XX 

Actually,  the  effect  of  the  pedestal  should  be  included  in  the  above  equations 
as  shown  in  Eq.  (12)  and  (13)  in  the  analysis. 
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The  numerical  method  uses  Eqs .  (1-a),  (2-a)  and  (3-a)  by  rewriting  them  into 
finite  difference  form: 


AV=  cu1($Aaz)-(x  +e) 


AM=  Vaz  +  <o*[(ir-ir)Az]  (&) 


Together  with  Eqs.  (5-a)  and  (6-a)  these  equations  form  a  set  of  recurrence 
relationships  which  can  be  solved  step  by  step,  starting  from  one  end  of  the 
rotor  until  reaching  the  other  end.  The  details  are  given  later. 


Occasionally  it  is  desired  to  perform  a  dimensionless 
governing  quantities  are: 


(7-a) 


(8-a) 


K-  = 


{Pi. 

(en. 


an a  lysis . 


The  two 


where: 


(EI)0-  reference  value  of  El,  lbs.in^ 

l  -  rotor  sPan  between  bearings,  inch 

My-  \fA4zj  total  rotor  mass,  lbs.sec^/in 

.  0 

Kc "  rotor  stiffness,  Ibs/in 

(<V  equal  to  or  proportional  to  a  critical  rotor  speed,  radians/sec 


For  a  uniform  shaft  (El  =  constant,  A  =  Constant): 
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where  fl  designates  the  order  of  the  critical  speed.  Thus,  for  the  first 
mode:  H  =  l  i.e.  , 


(ei)0=  -  ?.VH-er 


(Uniform  shaft,  first  mode) 


However,  it  is  not  necessary  that  be  a  critical  speed  but  Eq.(7-a)  must 
be  satisfied. 


The  dimensionless  parameters  become: 

x‘=  x/e. 

l/( 

(fi)’  =  er/ffi). 
v'  =  \/tX 
H1  *  m  /ari 

kjk 

UL)-<^lK- (^)(cw.)Civ') 
(jAV  -  t  fA/Mr 


where: 

e0-  reference  value  for  the  rotor  mass  eccentricity,  inch 
C,-  reference  value  for  the  radial  bearing  clearance,  inch 
C  -  actual  radial  bearing  clearance,  inch 
W#-  reference  value  for  tjie  bearing  reaction,  lbs . 

W*  actual  bearing  reaction,  lbs. 

L-  bearing  length,  inch 


The  dimensionless  bearing  coefficients  are  given  the  form  above  since  the  values 
obtained  from  lubrication  theory  are  CK^/VV  etc.  Normally,  a  di¬ 

mensionless  analysis  is  only  performed  for  a  simple  system  where  all  bearings 
are  identical,  i.e.  C~Cc  and  .  In  that  case  the  basic  dimensionless 

parameters  are: 

speed  ratio:  \<on/ 

dimen  ionless  ro^or  stiffness:  Xl~  Ckr/W 

dimensionless  bearing  coefficients:  C Kxir/W,  CujCm/Wj  (fc)  (CMty/wC)  ;eh. 
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Thus,  to  perform  a  dimensionless  calculation  for  a  given  value  of  kr  use  as 
input  to  the  computer  program: 

Speed  = 

.  hii  5 

Mass  at  station  l=  Jk  3.86069-10  (m=station  weight,  lbs; 

^  n=number  of  stations) 

(VU 


ftp-Ir) 


at  station  i  = 


Mr<a 


:±-  3.86069-105 


Cross-sectional  moment  of  inertia  for  section  i-a+l):  looo-I/Io 

Young  modulus  =  1 

Length  of  section  i-a+l)  =  Ufl 

Bearing  spring  coefficient  =  (SH 
Bearing  damping  coefficient?  ( ^^*) 


Unbalance  such  that 


:  tu„(oi.ln)  =  tuJo7.i*)r  £17 7.1 

i  •  7 


Then  the  computer  output  will  give: 

amplitude  =  ^  and  ^ 

bending  moment  =  M*  -  M/Co  M  r  M 

transmitted  force  =  (actual  force)/e0Kr  =  (actual  force) / ^ W k"/ 
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ANALYSIS  AND  DESCRIPTION  OF  THE  COMPUTER  PROGRAM 

The  remaining  part  of  the  report  describes  the  basic  analysis  and  gives 
the  detailed  instructions  for  using  the  computer  program: 

PN0011:  "Unbalance  Response  of  a  Rotor  in  Fluid  Film  Journal  Bearings"  for  the 
IBM  704  digital  computer.  This  program  calculates  the  rotor  deflect¬ 
ion  and  bending  moment,  the  pedestal  deflection  and  the  transmitted 
force  resulting  from  a  specified  rotor  unbalance.  It  differs  from 
conventional  programs  by  taking  into  account  the  variation  of  support 
flexibility  and  damping  along  the  whirl  path  of  the  rotor. 

The  supports  for  the  rotor  consist  of  a  fluid  film  bearing  on  a  ped¬ 
estal,  both  members  possessing  flexibility  and  damping  for  translatory 
and  rotational  motion.  The  flexibility  and  damping  are  linear  in  dis¬ 
placement  and  velocity  respectively,  the  proportionality  factors  den¬ 
oted  as  spring  and  damping  coefficients.  The  fluid  film  is  represent¬ 
ed  by  4  spring  coefficients  and  4  damping  coefficients  for  translatory 
motion  and  similarly  for  rotational  motion,  thus  allowing  for  coupling 
between  the  motion  in  two  mutually  perpendicular  directions.  The 
pedestal  has  no  such  coupling  and  is  represented  by  2  spring  and  2 
damping  coefficients  for  both  translatory  and  rotational  motion  with 
corresponding  pedestal  mass  and  mass  moment  of  inertia.  Hence,  each 
point  of  the  rotor  will  whirl  in  an  elliptic  path  around  its  steady 
state  position. 

In  audition,  the  program  includes  the  effect  of  gyroscopic  moment  and 
provides  for  couplings  in  the  rotor. 

THEORETICAL  ANALYSIS 


The  analysis  is  an  extension  of  the  Myklestad-Prohl  method,  see  Ref.  1, 
2  and  3.  The  rotor,  which  is  actually  a  continuous  system  with  an  in¬ 
finite  number  of  degrees  of  freedom,  is  replaced  by  a  finite  number  of 
lumped  masses  connected  by  weightless  springs.  The  computer  program 
calculates  the  vibrational  response  of  this  equivalent  system  exactly. 
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Thus  the  accuracy  of  the  results  depends  only  on  how  closely  the 
idealized  system  resembles  the  actual  rotor. 

Starting  from  the  left  end  of  the  rotor,  the  program  calculates  step 
by  step  the  bending  moment,  shear  force,  slope  and  deflection  along 
the  rotor.  Neglecting  the  shear  force  contribution  to  the  deflection, 
we  get  from  Fig.  1 : 

(1)  Mn*i  =  Mi  +  Ln  Vn 

9  ru I  =  9r\  +  CLn  Mn  +  bn  Vn 

^  Xrttt  =  Xn  +  Ln  On  +  Mu  +  dm  Vr\ 


The  program  assumes  El  constant  between  mass  points.  At  the  mass 
points,  the  forces  acting  on  the  rotor  are  introduced.  Four  contrib¬ 
utions  exist:  (1)  inertia  force,  (2)  unbalance  forces,  (3)  bearing 
reaction,  and  (4)  gyroscopic  moment.  In  general,  not  all  4  contrib¬ 
utions  apply  to  each  mass  point. 


Inertia  force.  The  rotor  performs  harmonic  vibrations  at  the  same 
frequency  as  the  rotational  speed.  Thus  the  inertia  force  is: 


AZX  7 

(8)  -  m  - —  =  7naJlX 

c)tl 

(9) 

Jiz  1 

Unbalance  forces.  To  allow  for  change  in  circumferential  position  of 
the  unbalance  along  the  rotor,  the  unbalance  is  given  two  components 
Utt  and  Uy  .  This  gives  rise  to  an  X  and  y  component  of  the  unbalance 
force : 

(10)  (v^-l4^i)unb=  OOzUx  COS  tat  -  u)ZUH  smart: 

(11)  (  V»f,n-  CoiuJb  +  U)Z  Ux  S  inu?t 

Bearing  reaction.  The  bearing  supports  have  flexibility  and  damping 
for  both  translatory  and  rotational  motion  of  the  rotor.  Since  the 
equations  for  the  two  types  of  motion  are  analogous,  only  the  equat¬ 
ions  for  translatory  motion  will  be  derived. 


The  bearing  support  is  shown  in  Fig.  3.  It  consists  of  a  pedestal 
with  mass  (  Mox^Moy  )  ,  supported  by  springs  (^x,)Ry)  and  dashpots 
(Ok,^  ).  There  is  no  coupling  between  the  X  and  y  direction,  i.e. 
no  transfer  impedance,  nor  between  the  translatory  and  rotational 
motion.  The  pedestal  supports  the  bearing  fluid  film  which  is  rep¬ 
resented  by  4  springs  and  4  damping  coefficients.  If  the  relative 
motion  between  the  journal  center  and  the  bearing  housing  is  denoted 
( X  )  y  ),  then  the  bearing  reaction  becomes: 

^  Cl'xn"  =  CxxX*  “  %  “  Cxy  ^ 

“  KyxX  ~  C(jxX  -*  Kyy  ^  ~  Cyy  y 

Setting:  . 

X  *  Xc  cos  U)t  +  X  s  s  t  w  orb 
y1  =  y  c-  cos  ut  4-  yj  sin  u)t 

we  get  from  Newton's  second  law  for  the  pedestal  mass: 
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(Kxx+  Ax"  CJlMox)Xc  +  b)(Cxx+6*)Xs  +K*j  L|c  +  U)Cxy(Js-(&x'UZKx)Xt  +  ^&xb 

-(j(cXXi^)x'c +kxxHr^x)x;-^^£/c  f  ^  ^---^xc+d-^M^Xs 

(13) 

+  ^xXk+CKl^+^-^AI^  Iji  - 

-iO&jxXc+tyx'rdCn+tf'j'Mfy+^tS^ 

Solving  the  equations  we  obtain: 

(v„-  (-AV&tXc  -4 VI, X, -4 V« ys)  c* u 

(14)  +  UVixXc-AV^Xs+ilVaxy 

('aVc^WiXs' A^^‘' cosl^ 

+  (  A  Vo(hXc-^Vc^5  ^  ^  "  A  S  s)  wt 

where : 

A  Vtxx  =  K*x  £  4-  tJ  C  k*  Cj  +  Kxij  <|  -4  Cxcj  r 

A  ybx  ®"  Kxx^  +(aJ  Cxx  f  *"  +  LJCkl j  ^ 

A  Vex  =  Kxxln  +  ^**L  +^*i  S  +  U)Cx^i 

(15)  AVdx  =-Kw  t  4-wCxxln  "Kxv,  t  +iJCx^S 

A  V<u^  =  K^x  K  4  UP  Cqx  l  +  klvfl  S  +  U)  Cy f  t" 

AVbij  -  -  Kvj*L+  t^Cijx  h  "  Kij'ft  +u>£^(j  s 

a  Vci^  =  f  4- 1 x;  C^x  ^  ^  +■  u)  Cij*|  ^ 

A  Vcit|  =-Kt^ 


and : 


r  GE+FD 

+  Fz+  62 
l  -  drJ+FH 

F  s  f\~  Kxy  Q  +  U)  Go  y  R 

H  -  -  Kx<j  S  +  co  C* j  T 

q  _  Q-  +  COGyvIn 

3  Fz+ 6* 

;  -  6H-FJ 

6  -  B  ~  R"(*> G<^  (^ 

3  *•  Kxy  T-  (O  Lk y  S 

ep  CO  CuxCL  -  Ku \*  b 

r\  2  1 1 » 

a  +  b 

r  <zc(  +  fee 

T>,C.  M  ZM, 

1*7 b1 

/4  *  Kxx  T  Ax  -0-)1MlX 

0*  a>  ('Cw  +(jj) 

0  *  $x -U)lMox 

E  *  to  6x 

0.  =  Kef y  +  -602  M/>y 

b  *  a;  (Ctfij  +  Gy  ) 

c/  * 

e  •  CO  Gy 

^  *  -  <Jf  -  Rg 

f  *  -  £  f 

s  -  S -Qh-  Rz 

t  s-T-tpA-t-R^ 

The  equations  for  rotational  motion  are  analogous  to  eq.  (14)except  for 
a  sign  reversal  (sign  convention,  see  Fig.  2): 

(*tAa*9c  +  AMb*6s  +AMcx<Pc  +  AMdx<Ps)  C05U)t 

+  i'AMbxdc  +  A  Max  d$  ~  AMix  <j Pc  ^  ^  Mcx  )  S\f\u>‘t 

(16) 

(Mi =  (*MCij9c  +AMd<j0s  May  Cfc  +  6  Miy  $s)  CoS  art 
+  9c  f  AMcy  9s  ~A  Mby  $  *  A  Mat/  )  Siniot 


where  the  coefficients  etc.  are  computed  from  eq.  (15) 

as  AtAur&Ya.*  ,  ^M(jx  —  etc.  by  replacing  the  translatory 

spring  and  damping  coefficients  by  the  corresponding  rotational  co¬ 
efficients. 


Since  the  fluid  film  coefficients  are  functions  of  speed,  directly 
through  the  Sommerfeld  number  and  indirectly  through  the  decrease  of 
eccentricity  ratio  with  increasing  speed,  the  computer  program  prov¬ 
ides  for  expressing  the  coefficients  as  a  function  of  speed,  e.g. 

4  ^  Kxx  =  Kxx,o  +  Kxx.i  *  U)  +  Kxx,z  *  40* 

and  similarly  for  the  other  coefficients,  (0  is  the  rotor  speed  in 
radians/ sec. 


Gyroscopic  Moment.  The  gyroscopic  moment  derives  from  the  change  of 
the  angular  momentum  vector  of  the  rotating  rotor  mass  as  it  whirls 
in  an  elliptical  path  around  the  steady  state  position  of  the  rotor. 
For  two  special  cases  the  gyroscopic  moment  is  known: 


(19) 


circular  whirl  path:  M6t(r.  -  (lp-IT  )  U) 

straight  line  (transverse  vibrations):  Mfiyr.  s  ~IT 


where  9  is  the  slope  of  the  rotor  deflection  and  Xp  and  It  are  the 
polar  and  transverse  mass  moment  of  inertia.  For  an  elliptical  path 
the  gyroscopic  moment  is  no  longer  linear  with  respect  to  the  slope  of 
the  rotor,  indicating  that  an  elliptical  path  is  actually  not  possible. 
However,  in  general  the  effect  of  the  gyroscopic  moment  is  not  too 
big  and  for  the  present  analysis  an  elliptical  path  will  be  assumed. 


The  coordinate  system  is  shown  in  Fig.  4,  where  0  is  the  steady  state 
shaft  center  position  and  0 '  is  the  whirling  shaft  center.  The  moving 
coordinate  system  )  is  defined  by  its  unit  vectors: 


e/ie'+v1  <#* 

Si +&*■+</>*  *  tfiTFr7ipL  *  tfTTerTifi- 


(20)  ei  “  si  </>*  i  J  el->-(pz  j  0 


'jJL=  jJL=, 

{(FTW1  r&  * 


-  —  £  ..  -—  £.  .  L  )  _  /  a  /«  j  \ 

»  {i+s^np-  j  ~  \  ’ r’  / 


The  angular  velocity  vector  becomes: 


•  i  #  «  .« 


<21)  <D-(«, ,«,.**)- 


The  moment  needed  to  sustain  the  motion  is  given  by  Eulers  equations1 

s  If  CO|  +  (  Zp~  IT)  Ll)i>  UJ,i 

(22)  =  ITW,  *  {I,-Tt)  W$0)j 

where  I  denotes  mass  moment  of  inertia  and  IT  and  If  -  Ip 

Let  us  first  assume  that  ( X  ?  \^  )  corresponds  to  the  directions  of  the 
major  and  minor  axis  in  the  elliptical  variation  of  the  rotor  slope. 


Therr: 


6,  *  £  cos  (uot+d.) 
(ft  -  G  5m  (ut  +  oc) 


Combining  eq.  (20),  (21)  and  (22): 


M,  -  -  It  «'«,  *  I,  ^ESjkK  J 


which  clearly  shows  that  the  gyroscopic  moment  is  not  linear  with 
respect  to  the  rotor  slope.  However,  only  the  first  harmonic  can  do 


work  on  the  rotor.  Hence  a  Fourier  analysis  will  be  performed.  The 
following  integrals  apply; 


<z  tt 


sinx  COSX  dx 
E*coslx  +6**1**/ 


=  0 


,itr 


sinx  coiX  olx 
(Btcos'lx  +62sir iV)1 


=  0 


s  )  V\  7x  dx 
E1cos1‘x+6i$in1X 


smV  dx 

(Etco$lx+6ish iV)1 


_ dX 

£lcos\*  +  6W/ 


_ gtt _ 

(Elcoslx+&liin7x)1 


—2JL 

&(E+ 6) 


_D_ 

Ed3 


zn 
r  6, 

rrfE'+t,1-) 

f36,3 


Then  the  first  harmonic  becomes: 


u‘<f, 

(24) 

My  =  ( Tp~ IT)  u>z9i 


In  the  limit,  eqs.  (24)  agree  with  eqs.  (19). 

Eqs.  (2  4)  must  be  transformed  back  to  the  actual  (0,<f) -coordinate 
system.  Setting 


&  =  9c.  Cos  U)t  t  sin  out 


(25) 


^  =  ^coos  Uit  +  SIM  K/t 


describing  an  elliptical  variation  of  slope,  we  get: 


I f  1/2  ♦«*)  * 


COS  2$  = 


l(o<l+6sH-tf)'+*M+tys)1 


.  ZfbA+Ml) 

p' 


where  (3  is  the  angle  from  the  position  X-axis  to  the  major  axis  E,  pos¬ 
ition  in  the  same  direction  as  U)  .  Then: 


6i  =  0  cos  ft  +  (f>siv\{i> 

(27) 

s*  -  (9  si*^  +  cj)  COS  jS> 

Substituting  eq.  (26) - (27)  into  eqs.  (24)  gives: 

(28)  M*  =  (Mxl  -  uoz[zaM6xIp~  &cIr] Cosu)t  Ho2[zAM^IpBsIT]s]^t 

(29)  -  Mx  =  ( My* =  ^t^^IP-fcrr]  CfiiuJt  f u/[2AM^rP-f  IT] sinu)fc 

where 

,,D.  ,,  (8c*ft)(0^s-«,C) 


(31) 


-  ^  rkj 

*'  (&♦*)*+ (fls-lfc)1 


Since  eq,  (28)  and  eq.  (29)  are  not  linear,  an  iterative  method  is 
used.  For  each  rotor  speed,  the  program  performs  a  number  of  iterat¬ 
ions.  The  first  iteration  is  done  without  gyroscopic  moment.  After 
the  first  iteration,  the  gyroscopic  moment  is  calculated  from  eq.  (28)- 
(29)  and  these  values  are  used  in  the  second  iteration  and  so  on.  The 
calculation  has  converged  when  the  relative  change  in  rotor  amplitude 
and  slope  between  two  iterations  is  smaller  than  a  specified  limit. 


EQUATIONS  FOR  ROTOR  CALCULATION 


The  bending  moment,  the  shear  force,  the  slope  and  the  deflection  are 
expressed  by: 

My  =  Mxc  co suit  +  Mxs  £>inu>t 
V*  "  Vxc  cos  oot  +  Vxs  smujt 

0  =6  c  co$  cot +■  6  s  siiA  tot 

X  —  Xc.  cos  u>t  +■  Xt 

and  similarly  for  the  \j -direction.  Then  eq.  (1),  (2),  (3),  (8),  (9), 

(10),  (11),  (14),  (16),  (17),  and  (29)  may  be  combined  to  give  the 
equations  used  in  the  rotor  calculation  (see  Fig.  2): 

Mxcn=  Mxcn  +4MAxn&n  +  4Ml>xh$Sft+  t  (Mxc*T  Mx^ayru 

M ysn  =  Mxsn~  &  Mbxn#cr,+  (Mxsn-Mxsi  )^o 

tAycy\~  ^  0c*\  Mav|*i  4!:u  4  -  VV 

Mys*  -  Mij ~ L  Mdtjn  +4  Mci^l 

V/tn  =  Vxc,rt-i"*’ [wAUJ  "dVax^J  Xcm"^  VfewXirt  “^VdXw^j $n  +  ul^Uxu 

(32)  Vx5w  ~  VxS,*-\  +  Xcn  +  VW*J  Vjh  ^4  \/jXV\i)cii\-2)\/£.Xn^iiA  “  tO  UlJia 

=  Vujc.n-1  -  Xc-vn  -  4  -^I4ynj  '  ii  y^rt  +C‘-)  tXijw 

yH,n  *  y^vi+^-n Xt*  -aVc^Xw  ^  4 -ojiuxa 

Mxc^fT  Mx£*i  +■  L*V*cw 

Myj,^i=  Mxs*  +  L*  Vx<a 


[Alj C.iwi"  M  ijOA  + 

0^**1*  Qcn  4  An  M*ca  +■  la* Vxlm. 

&it*.t~  diA  +  MtfS*  +  b^VyjA 

$:,nti  “  +■  6-a  b*Vw|4,* 

tfi  A+  ^  ^H'5*  4  UV4SK 
y^nti  *  Xtn  4  L.a0£*  4  Is*  Mxca  +  C^VycK 
Xj^-M  =  ^5A  4  L-U  0SA  +1oaMx»/\  "bC^Vx** 

”  Uc*  4  ^■v'  4  +0U  ^e* 

tfs,*»i=  tft*  4U$«+  UM^sa  4 

In  the  above  equations  An  j  6n ,  of*  are  given  by  eq.  (4),  (5)  and  (7), 

A^n - -dMdijn  and  aVo.m  ,  4/fcxKi - .iiUrfij*  by  eqs.  (15)  and 

(  Mx6*‘  Mxcm  )6tf;-- - )<|fW  by  eq.  (28)-(29). 

Eoundary  Conditions.  The  rotor  is  assumed  to  have  free  ends.  No  loss  in 
generality  occurs  by  this  condition  since  it  may  be  changed  by  letting 
the  end  points  have  bearing  support.  A  proper  choice  of  support  coef¬ 
ficients  will  then  allow  for  any  type  of  end  conditions. 

For  a  rotor  with  free  ends  the  bending  moment  and  the  shear  force  are 
zero  at  the  end: 


(33)  Myz.1  *  Mxn  =  Mtjti  s  Mijsi  ~  Vrc  1  ~  Vxsi 5  Vyc\  V/qji  =  0 

(34)  Mxcr  ~  lAxiv'  ~  Mijtr  *  M<jy  =  |/ycr  *  Vxsr  -  *  Xif*/-  =  0 
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Starting  from  the  left  end  of  the  rotor  (see  Fig.  2) ,  eq.  (33)  is  used. 
However,  the  slope  and  the  deflection  are  unknown.  Using  the  super¬ 
position  principle,  each  unknown  is  applied  separately.  A  summation 
gives  the  combined  effect.  Ten  calculations  are  performed,  using  eqs. 


(32) 

■ 

1. 

Ba  -  ! 

i =  (pa  ~  (flu  ~X'.i  *  x si  ~ ^ti* y i <  ■  o 

2. 

Bn  ~  l 

9ci  *  fa  —  Xci  —  X $ i  —  — 4jj,  =  Ui^—0 

3. 

<pc,"  1 

Qi.\  -  &<>  i  =fji  =  Xa  =  X*i  =  t|d|  =  t|S|  =  U  ^=Uq^=0 

4. 

<&«=  1 

=©si  =fci  =  Xti  —  X?i  =4cl”lx(ji  =Uyv\ 

5. 

Kcr  1 

0ci  =fo|  ~(fa  =  cfsi~Xsi  =  (|a  ~(>)}|  =  Wm\=1^«=0 

6. 

1 

&c\  ~Q\\  =  (flf i  =Xa  “^i— =^vk  —U^=0 

7. 

$6  1  ~  —  $  1  ~Xu  ~  Xjl~4$  1  =Mx*\-W‘ 

8. 

0o  =Qi—({cl  =$i — y.c\—  Xsi  =L/ti — 

9. 

—  Uxn  Lhfn 

&.  =  =*»=%=%  =  o 

10. 

Gyroscopic  moment  applied 

©ti  —  $$  i  =  i=  Xci  — Xj  i  —  (^i— (/<  i  =Ux*=Lh^rO 

For 

each  calculation  eqs.  (32) 

are  used  to  calculate  the  bending  mom- 

ent , 

the  shear  force,  the  slope 

and  the  deflection  along  the  rotor.  At 

the 

right  rotor  end,  station  r  , 

eq.  (34)  must  be  satisfied,  i.e. 

r 


(35) 


< 


fA/cr,i 

Mxcr^Z-- 

.  i 

- Micr.t 

—Mxis,  °\  -/VWj/fl 

tAnf.h- 

- Mviti 

— f  -fAti/,  \o 

i 

/ 

- M  t/cs,J 

(fin 

1  .  / 

tA^x— 

- Mi/y,  t 

>  < 

(fin 

>  =  < 

-  /vfi^(/o 

fari 

Zvur,  t  - 

- !urt  t 

N 

Yci 

~  Wcr,  <1  -  VV*./,  /  a 

- | /jfir  e 

Ym 

“  V^r,  f  -  l/ys  /,  i  e 

tw 

- | fy»t% 

y* 

■”U|6/(  q  ~  Vn  tv*,  i « 

IV.' 

- Xijir^g 

— 

*#« 

“Vijsr,  1  -l/q*', 

Eqs.  (35)  are  then  solved  for  Szi  - ffei,  and  the  actual  values  of 

bending  moment,  shear  force  etc,  along  the  rotor  can  be  determined. 

At  a  given  rotor  speed,  eqs.  (35)  are  first  solved  without  gyroscopic 
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moment,  i.e,  Mxcr^o" - ~  0  .  Then  the  gyroscopic  moment  is 

applied  according  to  eq.  (28) -(29)  and  new  values  are  found  for 

- Cfu  from  eqs,  (35).  This  process  is  repeated  until  at 

the  K'th  iteration; 


(36) 


ie!i- 


feiL  <fe.' 


0C-))| 


,..(*)  (Ml 

jiLHn  i!4i  I 


e“’l+l 


en+itfV - 1 a. 


(C) 

I 


where  is  the  convergence  limit  specified  by  the  computer  input. 

If  the  calculation  does  not  converge  within  a  specified  number  of  iter¬ 
ations,  the  program  goes  on  to  a  new  rotor  speed. 


In  the  computer  output,  the  rotor  deflection  is  given  by  the  dimensions 
of  the  elliptical  whirl  path.  We  have: 


X=XeC  OSUjt  +X«SlVwt 

(37) 

Lj=  Ljc  C.OSuO't  +  C|ssinU>t 

As  shown  in  Fig.  5,  the  ) -coordinate  system  is  rotated  an  angle  (3 

in  the  same  direction  as  u>  to  become  (Xi,Oi  ).  Then 


X,  =  Clcos  (cot ->-a) 

(38)  ' 

<j,  =  lo  +  <*■) 

where  CL  and  b  are  the  major  and  minor  axis  respectively  of  the  ellipse. 
F  rom  Fig.  5  : 


X,  -  XC.OS&  +  Cf 
Lji  =  -xsm/ 3  +Lj  C0Sj3 


Then: 

a. 

(39)  b 
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Here  it  is  necessary  to  allow  b  to  become  negative.  The  reason  is  that  the  ".rans- 
formation  from  the  x-y-coordinates  to  the  ellipse  must  be  able  to  discern  be¬ 
tween  forward  and  backward  whirl  (i.e.  the  shaft  center  may  travel  in  the  same, 
direction  or  in  the  opposite  direction  of  the  direction  of  rotation  depending 
on  the  values  of  x  ,  x  ,  y  and  y  ).  Let  the  angle  between  the  x-axis  ard 
the  instantaneous  radius  vector  be  y  : 


Then; 


y:  +ah*'(£) 


i.e. 


frcfc-XsVc)  >0 
(xcHs-Xs'lc)  <0 


•  {ofWArd  whirl 
'•  backward  wkii-l 

5  straight  liti;  ©i-brt  (b-o) 


Therefore; 

(39a) 


To  find  d  and  ^  expand  Eq.(38) 


(d) 

Then: 


(a)  ato5o<  =  Xc  loifi  4  lfc  Sihfi 

(b)  -a5ih<*  =  X}  (OS 6  4  Ifc  Sin/3 

(c)  bsrnd  =  ~X<5iti/3  -Hjccoiji 

bfoSoi 

c)  4  U)1 :  az+  b*  =  x*  4  Xjl  4  ij*  +  ijl 

(a)24(b)-fc ?-(<*)' :  al-b*  =  (xt+xrtf  + 

i.e.  cos 2/3=  sMfl* 

-U)-(\ ,)-  (cMd )  ;  i  fal-bv)  Sih2d  -  - (xcx5  +(/c^j) 

(a)1- (l)l4(c)’ -(</)':  (al-b7)  <os2«<  =  + 

i.e.  cos2o(  =  —4^"^ 

A  ~b  A  ~  b 

Thus  in  total; 


Next 


(40) 


toun  2  ft  - 


2(Xclt c  +  X<M^ 

Xc  tXt-lfS-Llt 


COS  E(S 


xl±xl^zHi 

V-V 


(41) 


to.*  Z<3  = 


Xc-Xs'^t,-Lj!’ 


cos  z<x  = 


xi-xt+HcW 

a>-W 


Thus  £  is  the  angle  from  the  positive  X-axis  to  the  major  axis  of  the 
ellipse  (positive  with^J  )  and  Ol  is  the  phase  angle  for  the  radius 
vector,  measured  positive  from  the  major  axis  in  the  direction  of  to  . 
The  computer  output  gives  O-fi)  {3  and  Ct  for  both  the  deflection  and  the 
bending  moment. 

Coupling  Stations.  The  programs  allow  for  couplings  in  the  rotor.  At 
these  stations,  the  bending  moment  vanishes,  i.e.  M\k0  (the  coupling 
point  is  taken  just  to  the  right  of  the  mass  station).  When  the  prog¬ 
ram  encounters  a  coupling  station,  say  station  L  ,  the  following 
equations  are  set  up: 


- M»d,i 

'ft,' 

Myteifi 

/  — 

- AAxcL'9 

Xci 

+Mxulm 

Mxji'i  - 

— 

©SI 

MxsLS 

—  Mxsii 

Xi  r 

_ 

4  M/sijio 

M <jd,i  X\^a  i  - 

— 

Cl 

MijdS 

- Afyct  * 

7“ 

Mydf  +  Myd'io 

—  /fysi.,*/- 

<pil 

or  upon  solving 

(43)  eL=aijX4+bL  (i.j-1,2,3,0 

where  9 1  *  9a,  9t  -  9 ci  etc.  and  X(  ■  etc.  Then  the  bending 

moment,  shear  force,  slope  and  deflection  before  the  coupling  station 
become  functions  of  Xc\)Xa)  ^je.i  and  only.  As  an  example,  let 
the  shear  force  at  a  station  be: 


Vxcr\  “  Vi  &ci  *  14  0si  +Vi<Pa  +Vs  Xa  +  VoXsi  +  Vi  ijci  +  1/f  +  lt0 


Introducing  eq. (43)  gives: 
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Xfcn  ^5  4  V,  4  (in  V2  4  ]  Vc,  4  [ V£ "+  V,  +  3*2  \^  4  C^jj  \j  4  ]  ^ 

+  +  4  ^33^3  +  tye,4  [  Vj  "*■  +  <?34  +^44^]^, 

4  4VI(+  \»,v,  +t=2v2+  k,V34 

and  similarly  for  V*Sfv  VJCh  ;VJSh>  MXCh/ - (jth. Instead  of  ©ci^i/fci  <W  fii 

we  have  as  new  variables  the  slopes  just  to  the  right  of  the  coupling 
station  Wy  l,C.  ©C(v /$*  1  (fi^heri  the  calculation  proceeds  as  be¬ 
fore  until  either  a  new  coupling  station  or  the  right  end  of  the  rotor 
is  reached. 

Transmitted  Force  and  Pedestal  Motion 

Let  the  force  transmitted  to  the  pedestal  at  station  n  be  denoted  F. 

From  Eq.  (12)  it  is  seen: 

^Vc  “  Mrc,h-/  ~  Mrcn  4  ^ ^Ch  4"  ^  ^xn 

(45)  FyS  -  VX5  h_t  “  Vx$h  4  CO  XJh  —  to 

Rjt  ~  ~  ^jcn  4  CO  LjCh  4-  £ol  Uijh 

f<j5  *  V,S|h.r  VV  +  "i-v'ijtH  •+  ^  ukh 

Denoting  the  amplitude  of  the  pedestal  masss  Xp  and  Ljp  (see  Fig.  3) 
we  get: 

X  Frt-iFyS 

p  ~  dex-MxtoM-ico<4x 

nD  =  — Kg  - 

j  if),  ~M^co*+ico^ 

or 

V  —  FXC  (Xk  ~  ^xtO*)  ~  Fys  CQ^y 
p£  (Xx-tW)Mcodx)1 

y  _  Fxc  cody  4  F»s  (Xy-Hxco1) 

?S  (Xx-Mx^^Ccody)1 

u  _  Fuc C ~  Co*)  —  Ftfj  Cod^ 

^PC  (^-M^Co^-fCcod,)1 

_  FweCodfc -+ FfsCiCd-MwCO1) 

^ 

The  force  transmitted  to  the  base  becomes: 

Py  =  *xXp  -HCOdy  Xf> 

l^j  '  Ljp  4  ICO  d’jj  tjp 

or : 


(47) 


fxc  =  4  *pC 

Pxs  ~  ^Vs  4  He  (*)  Xps 

~  r*/C  4  ^0  ^jpc 

R)S  r  R,S  4  Wl^s 

Energy  Balance 

Let  the  relative  amplitude  between  rotor  and  pedestal  mass  be  X/=X-Xp 
and  lj'~  at  a  bearing  station.  Then  the  energy  dissi¬ 

pated  in  the  bearing  and  the  pedestal  per  revolution  becomes: 

Energy  Dissipated  - 

(48)  1\{ cjCxx  Cxi'+xf)  4  (o»C^+co^)(xJ^+xJ  i/') 

“(k*V-  ky,)  (<  $  -  *J  4  tyi+Xf*  )  +  UJ^(ijpc+(jpS  )  } 
+'Jt{coDX!f(e,c2+a'1)  +c*>b^(fct+(f)sl)  +  ^Dxj4£o^*)(0<(P£+©j(pi) 

- (M«f  fy* ) (&$- e's  #)  4  *><** fejc  +efsh(jfy  (fit  +  fr\ ) } 

A  summation  over  all  bearings  gives  the  total  dissipated  energy. 

At.  each  unbalance  station  there  is  an  energy  input: 

(49)  Energy  Input:  7Tf  io%(xs-%)-KJUy(xs 

Summing  over  all  unbalance  stations  gives  the  total  energy  input 
which  must  equal  the  dissipated  energy. 

COMPUTER  INPUT 

The  input  data  is  prepared  according  to  the  following  instructions. 

Note  that,  unless  specifically  stated,  no  input  card  may  be  omitted. 

Card  1  and  2:  (72  cols.  Hollerith)  Identification:-  Any  descriptive 

text  may  be  punched  in  cols.  2-72.  These  two  cards  must  always  be 
included . 

Card  3:  (1015)  Control  parameters  - 

Word  1 .  Number  of  rotor  mass  stations  -  The  number  of  mass  stations  is 
determined  by  the  above  considerations.  Also,  there  must  be  a  mass 
station  at  each  rotor  end,  at  each  bearing,  at  each  unbalance  and  at 


each  coupling  point.  The  mass  at  a  station  may  be  zero.  The 
maximum  number  of  mass  stations  is  80. 

Word  2.  Number  of  bearings  -  This  integer  denotes  the  total  number  of 
bearings  along  the  rotor.  A  maximum  of  25  bearings  is  possible. 

Word  3.  Number  of  unbalance  stations  -  This  integer  gives  the  total 
number  of  mass  stations  at  which  unbalance  is  applied.  A  maximum  of 
80  unbalance  stations  is  possible. 

Word  4.  Number  of  coupling  stations  -  This  integer  gives  the  total 
number  of  coupling  points.  It  cannot  exceed  20. 

Word  5.  Pedestal  flexible/ rigid  -  If  this  integer  is  zero,  the  program 
assumes  the  pedestal  to  be  rigid  for  both  translatory  and  rotational 
motion  and  no  pedestal  data  is  included.  If  the  integer  is  1,  the 
pedestal  has  flexibility  and  damping  and  pedestal  data  must  be  provided. 

Word  6.  Support  tilting  -  If  this  integer  is  zero,  neither  the  bear¬ 
ings  nor  the  pedestals  resist  rotation.  In  that  case,  neither  the 
input  for  the  bearing  dynamic  coefficients  for  rotational  motion  nor 
the  pedestal  data  for  rotational  motion  can  be  included.  If  the 
integer  is  1,  the  bearings  and  the  pedestals  have  flexibility  and 
damping  for  rotational  motion. 

Word  7.  Gyroscopic  moment  -  If  this  integer  is  zero,  no  gyroscopic 
moment  is  included  in  the  calculation.  If  gyroscopic  moment  is 
desired,  the  integer  should  be  1. 

Word  8.  Number  of  computations  -  It  was  indicated  above  that  the  eight 
bearing  parameters  were  dynamic  coefficients  and  so  could  account  for 
the  variation  of  parameters  with  running  speed  in  an  approximate  manner. 
However,  if  a  more  precise  representation  of  these  parameters  is  des¬ 
ired,  these  values  can  be  entered  each  time  a  new  running  speed  is 
designated.  In  order  to  facilitate  this,  there  is  provision  in  the 


program  for  entering  only  the  bearing  or  bearing  and  pedestal  data 
and  the  corresponding  running  speed  without  re-entering  the  rotor, 
coupling  or  unbalance  data.  Then  this  word  8  of  the  control  para¬ 
meters  designates  the  number  of  sets  of  parameters  which  are  to  be  run. 

If  this  value  is  1,  the  program  assumes  that  the  bearing  data  is  being 
entered  as  coefficients  of  quadratic  equations  in  £0  .  Note  below  that 
the  input  format  of  the  bearing  data  differs  depending  on  whether  this 
value  is  equal  to  or  greater  than  one. 

Word  9.  Diagnostic  -  If  this  integer  is  zero,  no  diagnostic  will  be 

performed.  A  value  of  1  will  provide  the  diagnostic  output:  the 

diagnostic  increases  the  amount  of  output  a  considerable  amount  and 
is  provided  primarily  for  use  in  trouble-shooting  the  program  and  so 
this  value  should  always  be  zero. 

Word  10.  Input  -  If  this  integer  is  zero,  the  program  will  return  to 

read  in  a rgw  set  of  input  upon  completion  of  the  computation.  For 

the  last  set  of  input  this  value  should  be  1. 

Card  4.  (1P4E15.7) 

2 

Word  1  is  Young'  modulus  E  in  lbs/in  .  It  is  constant  throughout  the  rotor. 
Since  the  program  never  uses  E  by  itself  but  always  in  the  product  El 
(I -cross-sectional  moment  of  inertia)  any  actual  variation  in  E  can  be 
absorbed  by  changing  I  accordingly. 

Word  2  is  the  scale  factor  for  the  determinant  in  the  simultaneous  equation 
subroutine.  In  general  this  item  is  1.0.  It  is  a  factor  by  which  the 
determinant  is  multiplied  to  control  computer  over/under f low .  The  simul¬ 
taneous  equation  subroutine  is  used  4  places  in  the  program:  once  when 
solving  for  the  unknown  end  deflections  (i . e . Eq . (35) )  and  3  times  when  solv¬ 
ing  for  the  unkown  slopes  in  the  coupling  calculation  (i.e.  Eq.  (42)).  If 
an  over/underflow  occurs  during  the  calculation  the  program  output  will  con¬ 
tain:  "OVER /UNDERFLOW  IN  XSIMEQF  AT  _  _  (integer)"  where  the  value  of  the 
integer  is  1  to  4.  If  it  is  1,  2  or  3  the  error  is  in  the  coupling  cal¬ 
culation.  If  it  is  4  the  error  is  in  solving  Eq.  (35).  Changing  the  scale- 
factor  may  eliminate  the  trouble. 


If  the  determinant  is  singular  the  output  gives:  "MATRIX  IS  SINGULAR 

IN  XSIMEQF  AT  _ (integer)".  If  either  of  the  two  errors  occur  the 

program  proceeds  with  a  new  rotor  speed. 

ROTOR  DATA 

The  rotor  data  will  differ  depending  on  whether  the  effect  of  the 
gyroscopic  moment  is  included  in  the  computation.  For  the  case  where 
no  gyroscopic  moment  is  included;  i.e.  where  word  7  of  card  3  is  zero, 
the  rotor  data  is  entered  as  follows: 

Card:  (1P3E14.6)  -  An  input  card  must  be  given  for  each  mass  station. 

Each  card  has  3  items . 

Word  1  -  the  mass  at  the  station  in  lbs. 

Word  2  -  the  length  of  the  shaft  section  to  the  right  of  the  station 
in  inches . 

Word  3  -  the  cross-sectional  moment  of  inertia  of  the  shaft  section  to 
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the  right  of  the  station  in  in  . 

For  the  last  mass  station  the  shaft  length  and  the  cross-sectional 
moment  of  inertia  has  no  meaning  and  may  be  set  equal  to  zero. 

If  gyroscopic  motion  is  included  and  word  7,  card  3,  is  not  equal  to 
zero,  then  each  card  contains  two  more  items  in  addition  to  the  3 
items  indicated  just  above.  Also  for  this  case,  the  rotor  data  cards 
are  immediately  preceded  by  a  card  which  contains  two  values  defined 
as  follows: 

Card:  (15 , 1PE23 .6) .  Gyroscopic  moment  parameters  - 

Word  1  -  Number  of  iterations  -  For  each  rotor  speed  the  program  first 


calculates  the  unbalance  response  without  gyroscopic  moment.  Based  on 


the  thus  obtained  rotor  slopes,  the  gyroscopic  moment  is  computed  and 
applied  to  the  rotor,  resulting  in  new  values  of  the  slope  and  the  process 
is  repeated.  The  program  counts  the  number  of  iterations,  excluding  the 
calculation  without  gyroscopic  moment,  i If  the  count  exceeds  this  input 
item,  the  results  obtained  are  printed  out,  the  iteration  count  is  reset 
to  1  and  a  new  rotor  speed  calculation  starts. 


Word  2  -  Convergence  limit  -  After  each  gyroscopic  moment  iteration,  the 
following  relative  error  is  calculated; 

<r  i  ♦  i  9.y>i  - 1 1  & 

ie,i"“l+|s>ltoi  +  l<^)"l+— -t  l4sf“| 


where  Oc^&si  )$c.i  and  are  the  slopes  and  X<-< ,  ;  Wei  and  ^51  are 

the  deflections  at  the  left  rotor  end.  The  superscript  is  the  iteration 
number.  For  each  iteration  the  computer  output  gives  the  iteration  number 
and  the  error.  When  the  error  is  less  than  or  equal  to  the  input  conver¬ 
gence  limit,  the  program  prints  the  results,  resets  the  iteration  count  to 
1  and  proceeds  with  a  new  rotor  speed. 


Following  this  card  are  the  rotor  data  cards. 


Card :  (1P5E14.6).  An  input  card  is  required  for  each  mass  station.  Each 

card  contains  5  items;  the  first  3  words  are  the  same  as  those  for  the  non- 
gyroscopic  moment  case  above  and  the  remaining  two  are: 
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Word  4  -  the  polar  mass  moment  of  inertia  in  lbs. in 

2 

Word  5  -  the  transverse  mass  moment  of  inertia  in  lbs. in 
LOCATION  OF  BEARING  SUPPORTS 

Card :  (1415).  This  list  provides  the  numbers  of  the  mass  stations  at 

which  there  is  a  bearing. 


UNBALANCE  DATA 

Card :  (15,  1P2E15.7).  A  card  is  provided  for  each  of  the  unbalance 

stations.  Each  card  contains  3  values: 

Word  1  -  an  integer  which  denotes  the  number  of  the  mass  station  at  which 
the  unbalance  applies. 

Word  2  -  the  cosine  component  of  the  unbalance  in  oz.  in. 


Word  3  -  the  sine  component  of  the  unbalance  in  oz.  in. 


By  providing  two  unbalance  components,  it  is  possible  to  take  into 
consideration  the  circumferential  variation  of  unbalance  along  the 
rotor. 

COUPLING  DATA 


If  the  rotor  does  not  contain  couplings,  (Word  4,  card  3  is  zero),  then 
no  coupling  data  is  necessary.  If  word  4,  card  3  is  not  zero,  the 
following  card  must  be  included. 

Card:  (1415)  -  This  is  a  list  of  integers  denoting  the  mass  stations 

at  which  there  is  a  coupling. 

PEDESTAL  DATA 


If  the  pedestal  is  considered  to  be  infinitely  rigid,  then  no  pedestal 
data  is  required.  In  this  case  word  5,  card  3,  pedestal  flexible/rigid 
is  zero.  Otherwise,  pedestal  data  is  required.  The  pedestal  data, 
like  the  bearing  data,  is  separated  into  translatory  and  rotational 
parameters.  Also,  as  before,  the  control  parameter  is  the  word  6, 
card  3,  support  tilting. 

Card:  (1P6E12.4)  -  A  card  must  be  provided  for  each  bearing.  On  it 

are  6  values  as  follows: 

Word  1  -  the  weight  of  the  pedestal  in  the  X  coordinate  in  lbs. 

Word  2  -  the  pedestal  stiffness  along  the  X  coordinate  in  lbs/in. 

Word  3  -  the  pedestal  damping  along  the  X  coordinate  in  lbs-sec/in. 

Word  4  -  same  as  word  1  but  for  the  coordinate. 

Word  5  -  same  as  word  2  but  for  the  coordinate. 

Word  6  -  same  as  word  3  but  for  the  coordinate. 

If  word  6,  card  3,  support  tilting,  is  not  zero,  then  all  of  the  cards 

concerned  with  the  translatory  parameters  are  followed  by  the  cards  for 
the  rotational  parameters.  Again  there  are  6  values  on  each  card  as 
follows: 

Word  1  -  the  mass  moment  of  inertia  of  the  pedestal  mass,  associated 
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with  the  X  coordinate  in  lbs.in 


Word  2  -  the  pedestal  spring  coefficient  for  rotational  motion,  assoc¬ 
iated  with  the  X  coordinate  in  lbs-in/rad. 


Word  3  -  the  pedestal  damping  coefficient  for  rotational  motion,  assoc¬ 
iated  with  the  X  coordinate  in  lbs-in. -sec/rad. 

Word  4  -  same  as  word  1  but  associated  with  the  4  coordinate. 

Word  5  -  same  as  word  2  but  associated  with  the  4  coordinate. 

Word  6  -  same  as  word  3  but  associated  with  the  coordinate. 

SPEED  AND  BEARING  DATA 


Each  bearing  is  represented  by  16  dynamic  coefficients,  8  for  trans- 
latory  motion  and  8  for  rotational  motion.  Of  the  8  coefficients,  4 
are  spring  coefficients  and  4  are  damping  coefficients.  Since  the  co¬ 
efficients  in  general  change  with  speed,  each  coefficient  is  expressed 
by  three  components;  e.g. 

K™  *  K*x,o  +  W  +  Kxx,i  to2 

where  4)  is  the  rotor  speed  in  rad/sec,  Kxx,®  in  lbs/ in.  ,  Kxx.i  in 

2  2 

Ibs-sec/ in. rad.  and  in  lbs-sec  /in. rad  .  Similar  equations 

hold  for  the  other  15  coefficients.  As  indicated  earlier,  if  it  is 
desired  to  enter  the  bearing  data  at  each  value  of  frequency,  there  is 
provision  for  this  in  the  program. 

If  word  8,  card  3  is  1,  the  program  assumes  the  bearing  data  is  prov¬ 
ided  as  frequency  dependent  coefficients.  In  this  case,  a  card  is 
provided  with  the  speed  range  and  increment  and  this  is  followed  by 
the  bearing  data.  An  input  card  is  given  for  each  coefficient  at  each 
bearing.  Each  card  contains  three  items,  namely  the  above  mentioned 
three  speed  components.  The  sequence  of  the  input  cards  is:  first 
all  the  cards  for  the  translatory  motion  and  then  all  the  cards  for 
the  rotational  motion.  The  cards  for  the  rotational  motion  are  not 
required  if  word  6,  card  3,  support  tilting,  is  zero.  The  cards  should 
be  given  in  the  following  order:  |WfWj  < x*,  fyf  I****  fyx 

for  bearing  1,  Kvv  fU)CXx - caCyx  for  bearing  2,  etc.  to  the  last  bear¬ 

ing,  then  (if  word  6,  card  3  f  0)  ,  Mx* 


for  bearing  1,  etc.  to  and  including  the  last  bearing. 


Card:  (1P3E14. 6)  -  Speed  data. 

Word  1  -  initial  speed. 

Word  2  -  final  speed. 

Word  3  -  speed  increment. 

Card:  (1P3E14. 6)  -  Bearing  data  -  in  the  order  defined  above  with 

three  values  on  each  card  as  follows: 

Word  1  -  the  coefficient  A0  of  the  expression  /$*  /)0  +>9(  OJ  +  UJZ 
Word  2  -  the  coefficient  A(  of  this  expression. 

Word  3  -  the  coefficient  At  of  this  expression. 

If  word  8,  card  3,  is  greater  than  1,  the  program  assumes  the  bearing 
data  will  be  provided  for  each  value  of  speed.  In  this  case,  a  card 
is  provided  with  a  single  speed  value  and  this  is  followed  by  the 
bearing  data  as  follows:  all  of  the  translatory  stiffness  and  damping 
coefficients  are  provided  in  the  order;  two  cards  for  each  bearing. 

The  first  card  contains  the  X  coordinate  translatory  coefficients 
Kxx  }UiCyXj  and  Cu^and  the  second  card  the  coordinate  translatory 
coefficients  ky*  and  toCy*,  both  cards  for  bearing  one  followed 

by  two  cards  for  bearing  two,  etc. ,  to  the  total  number  of  bearings  . 
Again,  if  word  6,  card  3  is  zero,  no  rotational  parameters  are  required, 
otherwise,  they  are  provided  in  a  similar  manner:  one  card  of  values 
and  a  second  card  of  for  bearing 

1  followed  by  two  cards  each  for  the  remaining  bearings. 

The  card  format  in  this  case  is  then: 

Card:  (1PE14.6)  -  Speed. 

Card:  (1P4E14.6)  -  Bearing  data  in  the  order  defined  above  with  these 

values  on  each  card  as  follows: 

Card  1  -  Word  1  - 

Word  2  -coC 
Word  3  - 
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Card  2  -  Word  1  -  £44,^ 
Word  2  -cjC^,a 
Word  3  -  |CMX(A 
Word  U  -aC^^ 
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COMPUTER  OUTPUT 

The  computer  output  is  largely  self-explanatory  and  each  output  item  is  identified 
by  a  descriptive  text.  There  sample  cases  are  shown  on  Page  Alto  A17  .The  output 
first  lists  all  the  input  data,  i.e.  the  two  heading  cards,  the  control  words, 

Youngs  Modulus,  the  rotor  data,  the  bearing  stations,  the  unbalance  data,  the 
coupling  stations,  the  pedestal  data,  the  speed  data  and  finally  the  bearing 
data.  Thereafter  follow  the  results  of  the  calculations  with  one  set  of  output 
for  each  specified  rotor  speed.  First,  the  speed  is  given  in  RPM  which  may  be 
followed  by  the  intput  bearing  data  if  it  is  new  fpr  every  speed.  Next,  a 
9 -column  list  gives  the  rotor  amplitude  and  bending  moment  at  each  rotor  station. 
Both  the  amplitude  and  the  bending  moment  require  four  quantities  for  their 
description.  Since  each  rotor  station  whirls  in  an  elliptical  orbit  it  is  con¬ 
venient  to  express  the  four  quantities  in  terms  of  the  dimensions  of  the 
ellipse.  Then  the  four  quantities  become: 

1.  the  major  axis  of  the  ellipse:  a(i.e.  the  maximum  amplitude  or  the 
maximum  bending  moment  during  one  revolution. 

2.  the  minor  axis  of  the  ellipse:  b 

3.  the  angle  between  the  x-axis  of  the  overall  reference  system  and  the 
major  axis  of  the  ellipse:  (3,  degrees  (in  output  identified  by: 

ANGLE  X- MAJOR) 

4.  the  phase  angle  with  respect  to  the  cosine-component  of  the  unbalance: 

ot  ,  degrees. 

The  amplitude  is  given  in  thousands  of  an  inch  (mils)  and  the  bending  moment 
is  given  in  lbs. in. 

The  selected  method  of  presentation  is  illustrated  by  Fig.  5  and  is  given  in  de¬ 
tail  in  the  analysis  by  Eqs .  (37)  to  (41).  However,  a  general  description  will 
also  be  given  here. 

The  presentation  is  based  on  two  reference  coordinate  systems.  The  first  reference 
system  is  the  stationary  x-y  system  fixed  with  respect  to  ground,  and  which  has  at 
each  rotor  station  its  origin  in  the  center  of  the  statically  deflected  rotor 
(i.e.  the  deflection  due  to  gravity).  The  x-y-system  is  "couipunicated"  to  the 
rotor  via  the  specified  values  of  the  bearing  spring  and  damping  coefficients 
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(  KX)(  f  kyi^  f  UlC^ftc.  )  and  the  corresponding  pedestal  data.  In  other  words,  the 
directions  of  the  x-axis  and  the  y-axis  are  chosen  when  preparing  the  computer  in¬ 
put  and  the  choice  reflects  in  the  input  values  used  for  Ky*;  ete.  Then  the 
elliptical  rotor  orbit  is  centered  in  the  origin  of  the  x-y-system  (i.e.  the 
steady  state  shaft  center),  it  has  a  major  axis  a,  a  minor  axis  b,  and  the 
orientation  of  the  ellipse  is  defined  by  the  angle'  (3  between  the  x-axis  and 
the  major  axis,  measured  in  direction  of  rotor  rotation.  Note,  that  both  a,  b 
and  p  vary  along  the  rotor.  A  negative  value  for  b  signifies  backward  whirl. 

Thus  a,b  and  p  specify  the  dimensions  and  the  orientation  of  the  elliptical 
orbit  but  one  more  quantity  is  needed  to  specify  the  position  of  the  moving  shaft 
center  on  the  ellipse  at  any  given  time.  The  phase  angle  o(  serves  this  purpose. 
Let  the  major  axis  be  the  x^-axis  and  the  minor  axis  the  y^-axis  (see  Fig. 5), 
i.e.  the  x^-y^-system  is  obtained  by  rotating  the  x-y-system  an  angle  p  in 
the  direction  of  rotor  rotation.  Then  the  instantaneous  position  of  the  shaft 
center  is  given  by: 

y,  -  ac6$(tut+o<) 

Ljf  ~  b  Sih  (Cot  +ol) 

Note  that  the  orientation  of  the  x^-y^-system  changes  along  the  rotor  since  |3 
does  with  respect  to  the  x-y-system  the  instantaneous,  shaft  center  position  is 
given  by: 

X-'J  (acosfif+ibsih/!))1  cos  (lot*  *  + 

If  -  V  (b(oift)1  S'*  (cvt  + U  +  ta*  ft)) 

In  addition  to  determining  the  instantaneous  position  of  the  shaft  center  with 
respect  to  a  stationary  coordinate  system  it  also  may  be  desired  to  know  the 
position  with  respect  to  t.he  rotor  unbalance.  The  location  of  the  unbalance  in 
the  rotor  is  defined  by  a  coordinate  system  which  is  fixed  in  the  rotating  shaft 
and  whose  axes  are  called  "the  cosine  axis"  and  "the  sine  axis".  Hence,  each 
unbalance  consist  of  two  components:  a  cosine  component  and  a  sine  component 
(in  the  analysis  the  symbols  Ux  and  are  used,  respectively,  see  Eq.s(10) 
and  (11)).  The  instantaneous  orientation  of  the  cosine-s-ystem  is  defined  by  the 
angle  (  (oi  )  between  the  fixed  axis  and  the  cosine-axis.  Thus,  the  instantaneous 
phase  angle  between  the  amplitude  vector  (i.e.  the  radius  vector  from  the  center 
of  the  elliptical  orbit  going  through  the  instantaneous  shaft  center  position) 
and  the  total  rotor  unbalance  vector  is: 
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angle  by  which  amplitude  vector  lags  unbalance  vector  = 

'jb  ■itan'1  ( j  ~iah~'(a  iaii(tot+o()) 


Here,  £  ^Kh  and  indicates  the  summations  of  the  cosine-components 

and  the  sine-components,  respectively,  of  all  unbalances.  It  is  seen  that  the 
lag-angle  is  not  constant  as  the  shaft  center  moves  around  its  orbit,  It 
attains  its  maximum  and  minimum  values  when: 

Although  the  discussion  above  is  primarily  aimed  at  describing  the  motion  of 
the  shaft  center  (i.e.  the  computer  output  for  the  amplitude)  the  same  des¬ 
cription  applies  to  the  output  for  the  bending  moment.  However,  for  each  rotor 
station  the  output  lists  one  line  for  the  amplitude  but  two  lines  for  the  bend¬ 
ing  moment.  Whereas  the  output  for  the  amplitude  applies  at  the  rotor  station 
itself  the  bending  moment  has  one  value  immediately  to  the  left  of  the  station 
and  another  value  immediately  to  the  right  of  the  station.  The  output  gives 
the  left  hand  value  first  (i.e.  the  output  gives  Mn  and  M h  respectively,  see 
Fig.  2).  The  two  values  are  in  general  the  same  unless  the  particular  station  is 
a  bearing  station  which  resists  tilting.  The  last  listed  value  of  the  bending 
moment  should  always  be  zero  (i.e.  corresponding  values  of  the  major  and  minor 
axis  should  be  zero).  In  general  they  are  not  exactly  zero  but  very  small.  The 
amount  by  which  the  values  differ  from  zero  gives  an  indication  of  the  accuracy 
of  the  calculation.  Note,  that  for  this  reason  the  last  values  of  the  angles  f$ 
and  ei  are  meaningless . 
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Following  the  output  for  the  amplitude  and  for  the  bending  moment  come  the  re¬ 
sults  for  the  force  transmitted  to  the  bearing  housing  (equal  to  the  dynamic 
bearing  reaction).  If  the  pedestals  are  flexible,  the  force  transmitted  to 
the  foundation  and  the  amplitude  of  the  pedestal  mass  are  also  given.  Each  of 
the  three  quantities  are  presented  in  two  ways:  first  in  terms  of  the  corres¬ 
ponding  ellipse  (i.e.  in  analogy  to  the  rotor  amplitude)  and  secondly  by  their 
x  and  y-components .  Thus,  if  the  transmitted  force  is  F  the  output  gives  8 
quantities:  the  major  axis,  the  minor  axis,  the  orientation  angle  fit,  the  phase 

angle  ol  }  IFJ  j  |F^|  and  My  where  the  last  four  items  are  defined  by: 


force  in  x-direction:  Fx  -  IF*I  COS(b>t+o(x) 
force  in  y-direction: 


The  transmitted  force  is  given  in  lbs  and  the  pedestal  amplitude  in  thousands 
of  an  inch  (mils). 


The  next  line  of  output  serves  as  a  check  on  the  calculation.  It  gives  the 
energy  per  revolution  put  into  the  system  by  the  unbalance  forces  and  the  energy 
dissipated  per  revolution  in  the  bearings  and  pedestals.  Theoretically,  the 
two  values  should  be  equal  but  numerical  inaccuracies  in  the  calculations  always 
cause  discrepancy.  Normally  they  differ  in  the  fifth  or  sixth  decimal  place. 

The  energy  is  given  in  lbs . inch/revolution. 

To  convert  it  into  HP  multiply  the  energy  value  by  the  speed  in  RPM  and  divide 
by  3.96  •  105. 


If  the  input  does  not  include  any  gyroscopic  moment  effects  the  calculations  are 
repeated  for  a  new  rotor  speed  and  the  output  will  follow  the  description  given 
above.  If  the  gyroscopic  moment  is  included  there  are  two  sets  of  output  for 
each  rotor  speed,  each  set  having  the  format  as  explained  above.  The  first  set 
applies  to  a  rotor  without  any  gyroscopic  moments,  and  the  second  set  gives 
the  final  result  for  the  calculation  with  the  gyroscopic  moment  included.  The 
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two  sets  are  separated  by  a  two  column  list  giving  the  sequential  results  of  the 
iterations  needed  to  perform  the  gyroscopic  moment  calculation.  The  first 
column  identifies  the  iteration  number  and  the  second  column  gives  the  relative 
convergence  of  the  iteration  procedure  as  explained  in  describing  the  computer 
input . 
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FIG.  I  SHAFT  SECTION  BETWEEN  TWO  MASS  STATIONS 
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FIG.  2  CONVENTION  AND  NOMENCLATURE  FOR  ROTOR  CALCULATION 


FIG.  3  BEARING  AND  PEDESTAL  SYSTEM  FOR  TRANSLATORY  MOTION 


NOMENCLATURE 


A 

c 

Dm  &y*Ay 
E 


F,/v 


Ip 

IT 

KoCy  Kf^ 

Kr 

L 

L 

l 

M 

M.A 

mt 

K% 

P>A 


2 

Cross  sectional  area  of  shaft,  in 

Major  and  minor  axis  of  ellipse,  in  (or:  lbs,  lbs. in) 

Influence  coefficients  for  shaft  section,  see  Eqs . (4)  to  (7) 
Radial  bearing  clearance,  inch 

Bearing  damping  coefficient  for  translatory  whirl,  lbs. sec/in 

Bearing  damping  coefficients  for  conical  idiirl,  lbs . in. sec/radian 

2 

Youngs  modulus,  lbs. in 

Rotor  mass  eccentricity,  inch 

x-and  y-components  of  bearing  reaction,  lbs . 

Cross-sectional  moment  of  inertia  of  shaft, 

Ch  between  stations  n  and  (n+l)Jt  in^ 

2 

Polar  mass  moment_of  inertia  of  a  rotor  mass,  lbs. in. sec 
(in  input:  lbs. in  ) 

2 

Transverse  mass  moment  of  inertia  of  a  rotor  mass,  lbs. in. sec 
(in  input:  lbs.in^) 

Bearing  spring  coefficients  for  translatory  whirl,  lbs/in. 

Rotor  stiffness,  lbs/in 

Length  of  shaft  section  between  station  n  and  (n+1) ,  inch 
Bearing  length,  inch 
Rotor  length,  inch 

Bending  moment  (M  to  the  left,  M1  to  the  right  of  station  n) lbs . 

n  n  * 

2 

x  and  y-components  of  pedestal  mass,  lbs. sec  /in  (in  input:  lbs). 

2 

Total  rotor  mass,  lbs. sec  /in 

Bearing  spring  coefficients  for  conical  whirl,  lbs . in/radian 

2 

Mass  at  rotor  station  n,  lbs. sec  /in  (in  input:  lbs) 
x  and  y-components  for  force  transmitted  to  base,  lbs. 
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t 

V 

w 

z 

<* 

e/<p 

CO 

coh 


Time,  seconds 

2 

Cosine  and  sine-components  of  unbalance,  lbs. sec  (in  input: oz . in) 
Shear  force  (V  to  the  right  of  station  n);lbs. 

Bearing  reaction,  lbs. 

Components  of  the  rotor  amplitude,  inch  (in  output:  mils) 
Coordinate  along  length  of  rotor,  inch. 

Phase  angle  between  amplitude  radius  vector  and  unbalance 
see  Fig.  5,  radians 

Angle  between  major  axis  and  x-axis,  see  Fig.  5,  radians 

Pedestal  damping  coefficients  for  conical  whirl,  lbs . in. sec/rad. 

otx  du 

=  y  y  components  of  the  slope  of  the  deflected  rotor, rad. 

Pedestal  spring  coefficients,  lbs/in. 

Pedestal  damping  coefficients,  lbs. sec/in. 

Angular  speed  of  rotor,  radians/sec. 

Critical  rotor  speed,  radians/sec. 


Subscripts  and  Superscripts 

X  in  x-direction 

in  y-direction 

**y*WM  fi  rst  index  gives  force  direction, second  index  gives  amplitude 

direction. 

cosine  component 
sine  component 
applies  to  station  n 
pedestal 

relative  between  journal  and  pedestal 
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SAMPLE  CALCULATION  NO.l 


ROTOR  IN  RIGID  PEDESTALS  »NO  GYROSCOPIC  MOMENT 


3.1300000E+07  1.0000000E 


2  •  1 40000E  +  0 1 
3 « OOOOOOE  +  O 1 

1  •  25000GE+02 
1 .240Q00E+02 
3.820000E+02 

8  •  7  6  00CK3E+0  2 _ 

7  •  770 000 E +  0  2 

5 .  7COOOOE  +  Q2 
4.530000E+C2 
2 » 350000E+02 

3  •  000000E+02 
3.1 30000E  +  Q  2 
8.200000E+01 
3.200000E+01 
3.672000E+01 

2  .  82 1 OOOE+0 1 

9  •  043000E  +  0 1 
4.297000E+00 
4. 297000E+00 
7  •  852000E  +  0 1 
3.423000E+01 
9  •  820000E  +  0 1 


2.431000E+02 
3  «  300000E  +  0  2 
2.843000E+02 
6  •  370000E  +  0 1 
3  16  22 


4.09000QE+00 
5 . 3 1C00QE  +  Q0 
6  *  6  200  00E  +  00 
2.25UOOOE+GO 
5.910000E+00 
9 . 780000E  +  Q0 
7.190000 E+00 
8. 570000E+00 
7.190000  E  + 6  0 
5  . QOOOQQE+OQ 
6.6200  0  6  E  +  0  Q~ 
4. 870000E+00 
2.1 30000E+00 
3 . 1 90QOQE+QO 
7.420000E+00 
8 . 650000E+00 
1.237500  E  +  0 1 
JU.2  3J7  500E  +  00 
1.237500E+00 
_4_»  750000E+00 
5  •  320000E  +  00 
6.1  90000E  +  00 


6.5  10000E  +  00 
6.2  lOOOOE+QO 
4. 3 300 0  0 E  +  00 
0.  OOOOOoE+OO 

27--' 


0  0  1 


t-00 

)_  §  •  2  50000E  +  0 1 

3  4.880060E+61 

3  4.880000E+01 

j  1 .44Q000E+Q2 
1  3.02_0000E  +  02 

3  7.788000E+02' 

J  2  •  580000E+03 
1  8 . 380000E  +  02 

3  3.950000E+02 

3  3.850000E+02 

3  3j>_2 1QQ0QE+02 

T  ^6.130bd0E+b2 
1  1.930000E+02 

3“  6  •  800000E  +  0 1 
3  4. 180000E+Q1 

3  2  •  OOOOOOE  +  O 1 

3  5 .476Q00E+00 

f  3.1 20000E-01 
3  3 . 120000E-01 

3  3  •  120000E-01 

3  2.000000E+01 

3  T .G92000E+02 
3  8 . 540000E+02 


3  1 . 602000E  +  03 

3  8 . 360000E  +  02 

3  1 . 554Q00E+02 

3  O.OCOOOOE+OO 


1.0000000E+01  O.OOOOOOOE+OO 


1  . 

1  .463000E+06 


57100000E+03 
0  .  003000E+C0 
0  .  OOGGOOE+Od 
0. OOOOOOE  +  OO 
0 .OCOOOOE+OO 
O.OCOOOOE+OO 
0, GOCOOuE+OO 
0. OOOOOOE  +  OO 
TT.TTOJITOOE+CKr 
0  .OOOOOOE+OO 
0.OOOOWE  +  OO' 
0. OOOOOOE+OO 
0.  OOOOOOE  +  OO 
0. OOOOOOE+OO 
O.OCOOOOE+OO" 
0. OO^OOOE  +  OO 

o.ooqooo"E+qo 

0. OOOOOOE+OO 

o.ooooooe+ocT 

0. OOOOOOE  +  OO 
O.OOOOOOE  +  O'O 
0. OOOOOOE  +  OO 
0.00013  CUE + OTT 
0.  OOOOOOE  +  OO 


1 .000000E  +  03  5. 100000E  +  03  2.000000E+03 

1  •  542000E  +  06  O.OOOOOOE  +  CO  0. OOOOOOE  +  OO 
3 .36 2000E  +  06  0.00  0000 E  +  00  0.0 CO 00 0E  +  00 

-8.010000E  +  05  0. OOOOOOE  +  OO  0. OOOOOOE  +  OO 

1.302000E+06  O.OCOOOOE+OO"  0 . QOOOOOE+OO 
1 . 542000E  +  0  5  O.OCOOOOE  +  OO  0. OOOOOOE  +  OO 
1  .O4OO0OE  +  O6  0,  OOCOOuE  +  OO"  0  .'dOCdOOE  +  OO  " 
9  •  720000E  +  0  5  0. OOOOOOE  +  OO  0. OOOOOOE  +  OO 

IT 3(1213110^115 - OTWJUTTO E+W'  O.OOOOOOE+OO"" 

1.542000E+06  0. OOOOOOE+OO  O.OCOOOOE+CC 

37 3 62000E+06  Q.OOOOOOE  +  00  '  0 . OOOOCOE'+OO  ~ 
-8.010000E  +  05  0. OOOOOOE+OO  0. OOOOOOE  +  OO 

1 .302000E  +  06  0.  OCOOOOE  +  OO  0 . 0  O'OOOO'E  +  0  0 

1.542000E+05  0. OOOOOOE+OO  O.OCOOOOE+OO 

i  • 040000E  +  06  O.OCOOOOE+OO  "  0 . OOOOOOE  +  OO 
9.720000E+05  0 . OOOOOOE+OO  0. OOOOOOE+OO 

1 .302000E  +  06  0. OOOOOOE  +  OO  O.OOOOOOE+OO 

1  •  463000E  +  06  O.OOOOOOE+OO  O.OOOOOOE  +  OO 
3  •  296000E  +  06  0. OOGOOCE+OO  O.OOOOOOE+OO 

-1  .045000E  +  C6  O.OOOOOOE  +  OO  O.OCOOOOE  +  OO 
1 .781000E  +  06  O.OOOOOOE  +  O'O  "0  .  GO'OdOOE+UO" ' 
-1 . 197000E  +  05  O.OOOOOOE  +  OO  O.OOOOOOE  +  OO 
177 54WOr+0'6  O.OOoUOUE+015  O.GECrOOOF+dd  "" 
1.175000E+06  O.OOOOOOE+OO  0 .000000E+00__ 

+00 

1.463000E+06  O.OOOOOOE+OO  O.OOOOOOE+OO 
"  3.296  0TJT5F+ ff  6"“0n70  0T5FTJUF+i30  OTFCO'OFOF+TFO 
-1 .045000E  +  06  O.OOOOOOE  +  OO  O.OOOOOOE  +  OO 
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-1 .197000E+05  0 • 00G000E+Q0  0 .OOOOGOE+OO 

1.754000E+06  O.OOoOOOE+QQ  O.OOOOOOE+OO 
1.175000E  +  06  O.OOOOOOE  +  OO  "  U.OCOCOOE+OO 
1.781000E  +  06  0 .  OOOOOCE  +  OO  O.OOOOOOE  +  OO 
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SAMPLE  CALCULATION  NO. 2 

ROTOR  IN  FLEXIBLE  PEDESTALS  AND  GYROSCOPIC 


27  4  3  C 

3  •  1 300000E  +  07  l.C 
7  1.0C 

7.600000E+01  3.75 

_2_._1 400001+01  4.09 

3.000000E+01  5.31 

1.250000E+02  6.62 

1 . 24000QE+02  2.2*5 

3.820000E+02  5.91 

8 . 760000E+02  9.78 

7 . 7 70000E+02 _ 7.1 9 

5 .700000E+02  8.57 

4 . 5  30  000E  +  0  2  7.19 

2.350000E+02  '  5.00 
3.000000E+02  6.62 

3  •  1300001+02  4".  8  7 

8.200000E+01  2.13 

3.200 66  6  E + 0 1  3 Vi 9 

3  •672000E  +  01  7.42 

2 • 821 000E+0 1  8.65 

9.043000E+01  1.23 

4. 297000 E+ 00  “1.23 
4.297000E+00  1.23 

7 .852C00E  +  01  “  4.75 
3  •  423000E  +  0 1  5.32 

9  .82 OWOT+OY  6.19 
2.431 000E+02  6.51 

‘3. 3000001 +02  6.21 

2.843000E+02  4.33 

6 . 3  70 000 1+01  0.00 

3  16  22  27 

1  1.0000000E+U1 

10  1.0000000E+01 


0  1 _ 1_  _  1 

1  .OOOOOOOE  +  OO 
1 • OOOOOOE-03 
3 »  750000E+00  5.250 

4.090000E+00  4.880 

5. 3100 OOF +00 4 .880 
6.620000E+00  1.440 

2 • 2500G0E+Q0  3.020 

5.910000E+00  7.788 

9.7  80  0  0  0  E  +  00  2.580 

_7_._1 900001  +  00  8.380 

8.5700001+00  3.950 

7.1 90000E+00  3.850 

5  • OOOOOOE+OO  3.210 
6 .6  200  OOE  +  OO  6.130 
47 870000E+00  1.930 

2.1 30000E+00  6.800 

3.1 90000E+00  4.180 

7 .4200001+00  2.000 

8.650000E+00  “  5.476 
1.2  37500E+01  3.120 

1.2375001+00  3.120 

1  •  2  37  500E+00  3.120 

4. 7500 OOE+OO  2.000 
5.32U000E+00  1.092 

6. 190000E+0C  ‘8.540 

6.510000E+00  1.602 

6.2 10000E  +  00  ‘8.360 

4 • 3  30000E+00  1.554 

O.OOOOOOF+CO  0. 000 
27 

E+U 1  O.OOOOOOOE+OO 
E+0 1  1.0000000E+01 


MOMENT 
0  1 


5.2500001+01 
4.880000E+01 
4.8800001+01 
1 .4400001  +  02 
3  •  02000  0E  +  02 
7.788000E+02 
2 . 58000  0E  +  03 
8.380000E+02 
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INPUT  FOR  PN0011: 

UNBALANCE  RESPONSE  OF  A  FLEXIBLE  ROTOR  IN  FLEXIBLE.  DAMPED  BEARINGS 
Card  1  Text  Col.  2-72 

Card  2  Text  Col.  2-72 


Card  3  (1015) 


1.  NS.  Number  of  rotor  mass  stations  (  4s.  80) 

2.  NB.  Number  of  bearings  (  4s  25) 

3.  NU.  Number  of  unbalance  stations  (  —  80) 

4.  NC.  Number  of  coupling  stations  (  4s  20) 

■5.  NPST.  0:Rigid  Pedestal  1:  Flexible  pedestal 

6.  NMOM.  0:  No  bearing  resistance  to  moment  1:  Moment  resistance  in  brgs . 

7.  NGYR  0:  No  gyroscopic  moment  1:  Gyroscopic  moment  calculation 

8.  NCAL  1:  1st  type  of  bearing  data  input  ^2:  2nd  type  of  bearing  data 

9.  0:  no  diagnostic  1:  diagnostic  given  input. 

10.  0:  More  input  follows  1:  last  set  of  input 


Card  4  (1P4E15.7) 

2 

_  1.  E,  Youngs  modulus,  lbs /in 

_  2.  Scale  factor  in  simultaneous  equation  solution 


IF  NGYR  =  1 
Card  (15 , 1PE23 .6) 

_  1.  NIT.  Number  of  iterations  in  gyroscopic  mom. Calc. 

_  2.  Convergence  limit  for  gyroscopic  moment  calcul. 


ROTOR  DATA 


A- 13 


If  NGYR  =  0,  use  only  first  3  columns,  FORMAT  (1P3E14.6) 
If  NGYR  =  1,  use  all  5  columns,  FORMAT  (1P5E14.6) 

Give  one  card  for  each  rotor  station,  in  total  NS  cards 


Rotor 
Station 
(don' t 
punch) 


Station  Mass 


lbs . 


Length  of 
shaft  sec¬ 
tion 
inch 


Cross  sec¬ 
tional  Moment 

of  Inertia 
■!  „4 


Polar  Mass 
Moment  of 
Inertia 
lbs . in2 


Transverse  Mass 
Moment  of  Inertia 
lbs . in2 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 


Rotor  Stations  with  Rearing  Support 
(1415) 

Give  NB  items 


Unbalance  Data 
(15 , 1P2E15 . 7 ) 

Give  one  card  for  each  rotor  station  with  unbalance,  in  total  NU  cards 

_ .  1.  Rotor  station  number 

_ .  2.  X-component  of  unbalance,  oz.inch 

_ .  3.  Y-component  of  unbalance,  oz.inch 


A-19 


Rotor  Stations  with  Coupling 
(1415) 

Applies  only  if  NC  4  0.  Give  NC  items 


Pedestal  Data  for  Translatory  Motion 
(1P6E12.4) 

Applies  only  when  NPST  =  1,  Give  one  card  for  each  bearing,  in 
total  NB  cards 


Pedes .Mass 
x-direction 

lbs . 

Pedes .Stif fn. 
x-direction 

lbs/ in 

Pedes . Damping 
x-direction 

lbs .sec/ in 

Pedes .Mass . 
y-direction 

lbs . 

Pedes .Stif fn. 
y-direction 

lbs . / in 

Pede  s.  Damping 
y-di rection 

lbs . sec/ in 

Pedestal  Data  for  Tilting 
(1P6E12.4) 


Applies  only  when  NPST  =  1  and  NMOM  =  1.  Give  on  card  for  each  bearing, 
in  total  NB  cards. 


Mass  Mom. 
of  Inert, 
x-direction 
lbs . in^ 

Angular 

Stif fn. 
x-direction 
lbs . in/ rad 

Angular 

Damping 
x-direction 
lbs .in. sec/rad 

Mass  Mom. 
of  Inert, 
y-direction 
lbs . in2 

Angular 

Stif fn. 

y-*direction 

lbs.in/rad 

Angular 

Damping 
y-dire  ction 
lbs . in . sec/rad 

A-2  0 


Tvoe  1  Bearing  Data,  NCAL  =  1 
Speed  Data  (1P3E14.6) 

_ 1.  Initial  speed,  RPM 

_  2.  Final  speed,  RPM 

_  3.  Speed  increment,  RPM 

Bearing  Coefficients  for  Translatorv  Motion 
(1P3E14.6) 

Give  8  cards  per  bearing,  in  total  8*NB  cards.  Each  card  gives  one  co¬ 
efficient  in  the  form:  K**=ky*)e  +  l<’yi,/fCd  +  k'y*/2£ol  ,  Co + Cxxti  4J  +  Gx(2  Co* , 

- Kxv 

_  _  _  Co  Cyy 

_ _ kvy 

_  _  _  CoCy  y 

-  -  -  kyy 

_  _  ___________________  CoCyy 

_ ^y 

_  _  _  C oCyy 

_ K, 

_  _  _  Co  Cyy 

-  -  -  k*y 

-  -  -  CoCyy 

- -  -  Kyy 

_  _  _  CoCyy 

-  -  -  kyy 

-  -  -  CoCyy 
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Applies  only  when  NMOM  =  1.  Give  8  cards  per  bearing)  in  total  8*NB 
cards.  Each  card  gives  one  coefficient  in  the  form: 
fyrx  =  Mxx,o  +  Mxx,i  W  Mxxfi  to1  f  U)  Dxx  ~  ^x,o  **"  Dxx,i  U)  +  Dxx’,  1  to*  , 

_ M** 

-  -  Co  D*» 

_ 

-  CoDxij 

-  - - 

_  _  _  ^*0  Dljtj 

_ M,* 

_  _  _  <A* 
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Type  2  Bearing  Data,  NCAL  ~  2. 

Repeat  the  following  input  as  many  times  as  given  by  NCAL. 

Speed  Data  (1P4E14.6) 

_ .  Speed,  RPM 

Bearing  Coefficients  for  Translatory  Motion 
(1P4E14.6) 

Give  2  cards  per  bearing  with  4  coefficients  per  card,  in  total 
2*NB  Cards 

-  Kx*;CvC)(x;K)i^4;(*ij 

-  -  -  - K 

-  -  -  -  Kxif,<oG()r( 

-  -  . —  . .  - — - 


Bearing  Coefficients  for  Tilting 
(1P4E14.6) 

Applies  only  when  NM0M=1.  Give  2  cards  per  bearing,  in  total  2*NB  Cards 

-  -  -  - Me*,  ,  M<1J( 

-  -  -  - 

-  - -  -  - M^cvIUMx^ujD*,, 

- -  -  -  - 
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ABSTRACT 


Laboratory  tests  proved  the  noise  reduction  ability 
of  hydraulic -supports  when  applied  under  the  bearings  of  a 
high-speed  rotor.  The  results  were  compared  with  calculated 
values  as  well  as  with  tests  conducted  with  rigid-supports. 
Curves  presented  show  that  the  force  transmitted  to  the 
bedplate  was  reduced  by  a  factor  as  great  as  15  over  that  for 
rigid-supports . 
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INTRODUCTION 


An  experimental  investigation  was  conducted  into  the 
ability  of  a  hydraulically-supported,  pivoted-shoe,  journal  bearing 
to  reduce  structureborne  noise  resulting  from  rotor  unbalance. 


High-speed  rotating  machinery  onboard  submarines  represents 
a  principle  noise  source.  Because  of  the  demand  for  quiet  operation, 
equipment  such  as  main  propulsion  turbines  and  turbine-generator  sets 
are  built  to  stringent  manufacturing  specifications,  isolated  with 
resilient  mounts  and  balanced  using  the  most  refined  techniques. 

Still  a  relentless  search  for  more  silent  rotating  equipment  continues. 


As  part  of  a  "Study  of  Noise  Sources  in  800  KW  Turbine- 
Generator  Sets",  Hagamanl  analyzed  one  promising  concept  which 
involved  flexibly  supporting  a  rotor  directly  at  its  bearings,  close 
to  the  noise  source.  He  concluded  that  a  large  reduction  in  trans¬ 
mitted  force  appeared  attainable  with  sufficiently  soft  supports. 
Since  the  bearings  were  in  the  direct  transmission  path,  the  noise 
originating  in  the  rotor  must  pass  through  the  bearings  to  excite 
the  bedplate  and  surrounding  structure .  Therefore,  bearing  load 
variation  (transmitted  force)  has  a  direct  relation  to  the  structure- 
borne  noise. 
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Lund,  Sternlicht,  and  McHugh  *  demonstrated  with 
mathematical  analyses  and  laboratory  tests  that  the  oil-film  in 
hydrodynamic  bearings  had  spring-like  properties  and  thus  acted 
as  a  flexible  support  with  the  potential  of  reducing  the  noise 
transmission.  However,  drastic  reduction  seemed  possible  only  by 
using  an  external  support  system. 


Coil-springs  held  the  bearing  in  a  study  conducted  by 
Fistere  and  Dickson^  to  examine  the  feasibility  of  the  flexibly- 
supported  bearing  for  submarine  propulsion  turbines.  Since  ship¬ 
board  operation  required  correcting  for  the  variation  of  bearing 
load,  which  would  cause  excessive  deflections  because  of  the  low 
spring-value  needed  for  isolating,  accent  was  placed  on  the  develop¬ 
ment  of  an  automatic  positioning  system.  Data  gathered  from 
extensive  laboratory  testing  proved  the  soft-support  scheme 
practical  and  an  excellent  method  of  reducing  structureborne 
vibrations . 


1 


Numbers  placed  superior  to  the  line  refer  to  the  Bibliography. 
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The  present  contract  advanced  the  knowledge  of  the  flexible- 
support  concept.  The  analytical  work  was  explained  in  other  parts 
of  this  report.  This  portion  is  devoted  to  the  experimental  investi¬ 
gation,  which  had  two  principle  purposes: 

1.  To  determine  the  feasibility  and  desirability  of  hydraulic- 
supports  for  isolating  a  rotor-bearing  system. 

2.  To  substantiate  the  analytical  tools  developed  under  this 
contract . 

Compared  to  mechanical  springs,  hydraulic -supports  are 
potentially  space-savers  since  accumulators,  which  can  be  placed 
remote  from  the  bearings,  provide  the  flexibility.  As  an  additional 
advantage,  the  servo-return  system  can  be  integrated  into  the 
hydraulic-supports . 

The  laboratory  model  was  not  designed  as  a  shipboard 
prototype,  but  emphasis  was  placed  on  determining  valid  and  mean¬ 
ingful  measurements  of  the  isolation  ability  of  the  supports.  No 
servo-return  system  was  included  since  this  was  examined  in  the 
previous  report^. 

With  a  known  unbalance,  an  unbladed  SS(N)585  propulsion 
rotor  was  spun  to  5000  rpm  on  both  hydraulic-supports  and  rigid- 
supports.  Curves  in  this  report  show  transmitted  force  levels  for 
each  configuration  and  comparisons  are  made  with  calculated  values . 

In  order  to  complete  the  dynamics  picture,  plots  giving  the  journal 
vibration  are  also  included. 

At  4500  rpm,  transmitted  force  (measured  in  lbs.)  dropped 
by  a  factor  of  15  on  hydraulic -supports  when  compared  with  rigid- 
supports.  However,  the  level  failed  to  go  as  low  as  predicted 
because  of  friction  in  the  hydraulic- supports . 

In  contrast  to  the  hydraulic-supports,  a  system  incorpor¬ 
ating  mechanical  springs  (Ref .4)  could  be  designed  to  overcome  the 
friction  disadvantage.  Because  of  the  establishment  of  profound 
noise  reduction  by  a  flexibly-supported  bearing  in  this  investi¬ 
gation  and  Ref.  4,  it  is  recommended  that  a  prototype  be  designed 
and  applied  to  submarine  high-speed  machinery,  but  with  mechanical 
springs  rather  than  hydraulic-supports. 
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TEST  APPARATUS 

To  determine  the  advantage  of  hydraulic -supports,  two 
test  build-ups  were  investigated:  (1)  rigid-supports  and  (2) 
hydraulic -supports . 

Assembly  drawings  are  included  in  Figs.  2  and  3,  while 
Fig.  1  pictures  the  arrangement  of  the  test  components.  A  small 
steam  turbine  powered  the  test  rotor  up  to  5000  rpm.  The  rotor, 
better  viewed  in  Fig.  4,  came  from  a  Skipjack  Class  SS(N)585 
submarine  propulsion  turbine. 

It  is  seen  in  Fig.  4  that  care  was  taken  not  to  interfere 
with  the  isolating  function 
of  the  hydraulic- supports . 

For  example,  flexible  hosing 
supplies  and  drains  lubri¬ 
cating  oil;  a  30  in.  long 
1-1/4  in.  diameter  jackshaft 
connects  the  drive  turbine; 
and  a  35  lb/in.  spring  (see 
sketch)  allows  adjustment  for 
bearing  torque. 

Since  pivoted-shoe 
bearings  are  inherently  stable, 
they  were  utilized  to  preclude 
the  possibility  of  oil-whip  due 
to  the  soft  supports.  A  draw¬ 
ing  of  the  shoes  is  reproduced 
in  Fig.  5 ?  while  Fig.  6  is  a 
photograph  of  a.  front  and  back 
view.  Each  60°  sector  was  supported  on  a  spherical  button  so  as 
to  have  rotational  freedom  in  both  the  circumferential  and  axial 
directions.  The  axial  feature  permitted  good  alignment. 

Instead  of  the  equal  spacing  of  three  or  more  shoes  around 
the  circumference  as  in  conventional  practice,  only  two  shoes 
cradled  the  rotor  as  shown  in  Fig.  7.  This  had  the  advantage  of 
eliminating  the  intermediate  preload  effects  of  the  upper  pads,  and 
simplified  the  force  transmission  path  for  a  more  accurate  analytical 
description.  Navy  2190TP  oil  (MIL-L-17331)  lubricated  the  shoes, 
with  the  oil  supply  temperature  maintained  within  t  2°F  with 
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Bailey  Meter  regulating  equipment. 
The  oil  flooded  the  bearing 
cavity,  with  the  flow  or if iced 
through  the  clearnace  space  of 
the  end  seals  as  illustrated 
in  the  sketch. 
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A  schematic  of  the 
hydraulic  system  is  shown  in  the 
adjacent  diagram.  Two  such 
systems  supported  each  bearing. 

The  bearing  housing  was  floated 
on  the  actuator  pistons,  with 
accumulators  included  in  the 
hydraulic  lines  to  provide  the 
flexibility.  Because  of  its  low 
viscosity,  kerosene  was  used  as 
the  hydraulic  fluid  to  reduce 
damping  effects .  A  Spraque 
diaphragm  pump  pressurized  each 
system  separately.  The  valve 
to  the  pump  could  be  shut  form” 
ing  a  closed  system  between 
the  actuator  and  the  Greer 
blatter-type  accumulator . 

The  hydraulic 

actuators  held  the  bearing  hous¬ 
ing  as  described  in  Fig.  8. 

Since  the  closed  system  appli¬ 
cation  demanded  zero  leakage,  a 
rolling  diaphragm  sealed  the 
piston  of  the  actuator,  an  RKL 
A14-0081  Rollomotor .  The  supports 
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were  located  along  mutually  perpendicular  axes  and  rested  on  Endevco 
force  gages.  The  force  transmitted  to  the  bedplate  was  completely 
described  by  the  transducers .  The  bearing  housing  sat  on  needle 
bearings  to  eliminate  tangential  forces .  In  order  to  reduce  the 
friction  to  a  very  low  level,  provision  was  made  to  rotate  the  inner 
race  of  the  antifriction  bearing  with  gear-motors.  However,  this 
refinement  was  unnecessary  when  the  transmitted  force  proved  larger 
than  that  predicted. 

Fig.  9  provides  a  close-up  view  of  the  hydraulic -support . 
The  safety  supports  (shown  hiding  the  contact  between  the  needle 
bearing  and  the  bearing  housing)  held  the  housing  when  the  hydraulic- 
supports  were  inoperative.  When  the  hydraulic  system  was  pressurized, 
a  gap  of  .040  in.  was  maintained  between  the  safety  member  and  the 
bearing  housing. 

In  order  to  provide  a  standard  for  comparison,  rigid- 
supports  replaced  hydraulic -supports  in  the  second  test  configura¬ 
tion.  The  plate  supports,  seen  in  Fig.  10  and  11,  simulated  rigid- 
supports  and  yet  allowed  the  force  to  be  separated  into  X  and  Y 
components  for  meaningful  measurements . 
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INSTRUMENTATION 


Table  I  lists  the  instrumentation  used  during  the  tests 
and  Fig,  12  shows  the  instrumentation  arrangement. 

Transducers,  to  monitor  transmitted  force  and  journal 
vibration,  were  concentrated  at  the  bearing  opposite  the  driven 
end  of  the  test  rotor.  Both  measurements  were  made  along  mutually 
perpendicular  axes,  45°  from  the  vertical. 

Model  2106  Endevco  force  gages ,  with  a  Model  2622  Endevco 
power  supply  and  Model  2616B  Endevco  amplifiers,  provided  a  voltage 
proportional  to  the  transmitted  force.  The  wave  shape  and  level 
were  determined  with  a  Hewlet  Packard  oscilloscope.  Model  122A. 

The  accuracy  of  the  oscilloscope  was  checked  with  a  calibrated 
Model  400D  Hewlet  Packard  voltmeter  and  found  to  be  within  the 
least  count  of  its  scale.  The  gages,  which  incorporated  a  crystal 
material  to  measure  the  dynamic  force,  had  a  stiffness  of  at  least 
2  x  l(r  lb/in.  The  sensitivities,  which  had  been  factory  determined, 
were  corrected  for  system  capacitance  and  the  mounting  bolts. 


Vibration  was  measured  with  Westinghouse  HQ  balancing 
instrumentation  composed  of  (1)  pick-up,  (2)  sine-wave  generator, 
(3)  amplifier,  and  (4)  meter  box.  The  sine-wave  generator  attached 
to  the  end  of  the  test  rotor  providing  a  signal  at  the  rotational 
frequency  while  pick-ups,  which  were  in  contact  with  the  journal, 
produced  an  output  proportional  to  velocity.  The  amplifier  enabled 
a  1.X,  3X,  and  1QX  amplification  of  this  signal.  The  pick-up  and 
generator  signals  were  multiplied  in  the  meter  box,  acting  as  a 
wattmeter,  to  provide  a  read-out  of  the  vibration  level.  By  this 
means,  only  the  operating  speed  component  of  vibration  was  dis¬ 
played,  and  higher  harmonic  were  eliminated  since, 
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The  vibration  equipment  was  calibrated  as  a  system  at  various  speeds 
using  an  eccentric  cam. 


An  Electro-Products  magnetic  pick-up  sensed  speed  from 
a  60  tooth  gear  and  an  E-put  meter  displayed  the  speed. 
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Spring  contact,  iron-constantan 
thermocouples.  Fig.  5,  sensed  the 
babbit  surface  temperatures  in  the 
bearing  pads.  These  temperatures 
were  measured  using  a  bridge- type 
Leeds  Northrup  potentiometer  with 
the  set-up  shown  on  this  page.  All 
thermocouples  had  been  calibrated  in 
temperature  baths . 

Fischer  Porter  flowraters 
indicated  the  oil  flow  to  each 
bearing. 
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TEST  PROCEDURE 


For  the  two  test  configurations:  (1)  rigid -supports  and 
(2)  hydraulic -supports,  transmitted  force  and  journal  vibration 
data  were  collected  for  a  known  unbalance  in  the  balance  plane. 

The  same  unbalance  was  inserted  180°  from  its  original  location, 
and  the  test  was  repeated.*  In  addition,  test  data  were  taken  with 
all  inserted  unbalance  removed,  leaving  only  residual  unbalance. 
The  values  were  combined,  as  described  in  the  Calculation  Section, 
to  eliminate  the  effect  of  the  residual.  Two  replicants  were  run 
of  every  test  to  guarantee  repeatable  results. 


Before  testing  it  was  intended  to  balance  the  rotor; 
however,  the  residual  unbalance  proved  sufficiently  low  to  forgo 
balancing. 


7 !/r6/'r?C. 


Sixteen,  equally  spaced,  1  inch  threaded  holes  composed 
the  balance  circle,  whose  location  is  noted  in  the  above  drawing. 
Each  balance  plug  weighed  approximately  160  gra. 


*-This  refinement  was  dispensed  with  in  the  case  of  the  hydraulic- 
supports  because  the  residual  unbalance  was  small  compared  with 
the  accuracy  of  the  experiment. 
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The  table  below  lists  the  nominal  conditions  of  t.he  test" 
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The  signals  from  the  two 
Endevco  force  gages,  located  at 
right  angles  as  shown  in  Fig, 8, 
were  displayed  on  the  corres¬ 
ponding  X  and  Y-coordinates  of 
the  two  channel  oscilloscope. 
The  peak-to-peak  values  of  the 
transmitted  force  (see  sketch) 
were  recorded,  although  single 
rather  than  double  amplitudes 
are  presented  in  the  curves  to 
be  discussed  later. 

Westinghouse  HQ  instru¬ 
mentation  measured  journal 
vibration  along  the  same  axes 
as  the  force  transducers . 

Again  peak-to~peak  values  were 
recorded. 


So  as  to  maintain  a  consistent  variation  of  the  Sommerfeld 
number  with  speed,  the  bearing  oil  flow,  supply  temperature,  drain 
temperature,  and  bearing  housing  pressure  were  monitored.  Measur¬ 
ing  the  bearing  pad  temperature  enabled  calculation  of  an  average 
viscosity. 


Some  measurements  were  characteristic  only  to  the  hydraulic- 
supports.  For  example,  the  precharge  air  pressures  in  the  blatter - 
type  accumulators  were  checked  before  and  after  each  day's  testing. 

If  this  value  varied  from  the  predetermined  value,  the  accumulator 
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was  recharged  before  further  testing.  In  addition,  the  hydraulic 
line  pressure  was  recorded  at  each  speed  point. 


Before  the  rotating  tests, 
the  stiffness  of  each  hydraulic- 
support  was  evaluated  using  dead 
weights  (see  sketch) .  The 
weights  were  increased  to  160  lb. 
in  20  lb.  increments.  A  dial 
indicator  in  line  with  the  pull 
measured  the  bearing  housing 
displacement.  The  pressure 
changes  in  the  hydraulic  lines 
were  also  recorded. 
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RESULTS 

1 .  Rigid-Supports 

Figs.  14  and  15  plot  the 
experimental  values  of  journal 
vibration  and  transmitted  force 
produced  by  a  unit  unbalance  at  the 
balance  plane.  In  Vol.  1  of  this 
report,  it  was  shown  that  the  oil 
film  properties  are  equal  in  the 
X  and  Y-directions  and  the  force  and 
vibration  responses  would  theoreti¬ 
cally  be  the  same  for  the  two 
directions.  Therefore,  the  X  and  Y 
components  of  the  test  values  were  averaged  for  plotting  in  Figs.  14 
and  15.  Single  amplitude  values  are  given  (see  sketch). 

For  reference,  the  unattenuated 
force  is  included  in  Fig.  15. 

The  unattenuated  force  is  defined 
as  the  reaction  at  the  bearing 
if  the  unbalance  is  applied  as  a 
non-rotating  force.  The  experi¬ 
mental  points  fell  close  to  the 
unattenuated  curve  although  both 
the  bearing  oil- film  and  rotor 
have  some  flexibility.  Resonance 
conditions  appeared  at  approxi¬ 
mately  2400  and  4200  rpm. 

Table  II  summarizes  the  test  data  for  each  run.  The 
calculated  values  of  Sommer f eld  number  for  both  hydraulic  and  rigid- 
supports  are  presented  in  Fig.  13. 

2 .  Hydraulic -Supports 

For  ease  of  comparison,  the  vibration  and  force  curves  are 
plotted  in  Figs.  14  and  15,  together  with  the  rigid-support  curves . 
Although  the  journal  vibration  amplitudes  are  near  the  same  values 
for  the  two  cases,  the  transmitted  force  levels  were  appreciably  less 
using  hydraulic -supports .  For  the  higher  speeds,  the  force  values 
level  off  and  remain  a  constant  magnitude  while  for  the  rigid- 
supports,  the  force  continues  to  increase  with  speed. 
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In  Fig.  16,  the  advantage  of  the  hydraulic-supports  is 
seen  by  plotting  speed  vs  the  force  ratio.  The  force  ratio  is 
defined  as  the  force  transmitted  with  rigid -supports  divided  by  that 
with  hydraulic-supports.  The  most  improvement  was  obtained 
where  the  force  ratio  was  15 . 

Table  III  presents  a  complete  summary  of  the  test  data  for 
hydraulic -supports . 


Difficulties  with  friction  in  the  hydraulic  system  appeared 
when  the  spring  value  was  checked  with  dead  weights,  as  described 
in  the  Test  Procedure  Section.  The  motion  of  the  bearing  housing 
was  less  than  expected  for  the  range  of  weights  applied. 


The  spring  and  piston  type 
accumulators*  originally  used 
were  found  to  be  the  chief 
sources  of  friction.  A  rolling 
diaphragm  or  Bellofram  sealed 
the  piston  for  zero  leakage. 

As  weights  were  applied,  the 
pressure  gage  responded  (see 
sketch),  although  the  dial 
indicator  registered  a  large 
degree  of  sticking  within  the 
accumulator „ 

Greer  blatter -type  accumulators 
were  substituted  because  of 
their  inherent  freedom  from 
friction  other  than  the  small 
amount  due  to  the  elastic 
extension  of  the  blatter. 


By  placing  weights  on  the  bearing  housing  and  noting  its 
displacement,  the  remaining  friction  was  evaluated.  The  pressure 
in  the  hydraulic  system  was  recorded  so  that  the  force  change  on 
the  actuator  piston  could  be  compared  to  the  weight  added.  As 
illustrated  in  Fig.  17,  which  presents  typical  results  for  one  of 
the  hydraulic-supports,  a  load  of  160  lb.  produced  a  force  of  116  lb. 


*-These  were  the  same  accumulators  described  in  detail  in 
Reference  4. 
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on  che  piston.  This  indicates  static  friction  between  mating  parts 
preventing  relative  motion  that  might  have  occurred  in  the  anti¬ 
friction  bearings,  grease  lubricated  bushings,  actuator  guides,  or 
the  actuator  rolling-diaphragm.  Fig.  18,  however,  indicates  little 
remaining  stiction  in  the  accumulator.  The  spring  constant  of 
4730  lb/in  deduced  from  this  curve  agrees  favorably  with  the 
calculated  value  of  5560  lb/in. 
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CALCULATIONS 


1.  Bearing  Load 

The  bearing  reactions  R  were  determined  from  the  rotor 
section  properties  given  in  Table  IV  by  assuming  that  the 
bearings  carried  the  weight  of  sections  1  to  17. 

The  bearing  diameters  D  and  the  lengths  L  listed  are 
measured  values. 


Turbine  End 
.19 

L  »  ~*32  *  3.594  in. 

D  -  6.9837  in. 

R  =  2161.2  lb. 

P  =  “  86.11 

LD  In^ 


L 

D 

R 

P 


Free  End 


3j|  -  3.594  in. 

6.9844  in. 


2293.2  lb. 


91.36  lb. 

5* 


2.  Bearing  Clearance 


The  clearance  was  found  from  the  difference  in  micro¬ 
meter  measurements  of  the  journal  and  bearing  shoe  diameters. 


Turbine  End: 


6.9978  -  6,9837 
2 


.00705  in. 


Free  End: 


6.9978  -  6.9844  .  .0067  in. 


3.  Viscosity 


The  viscosity-temperature 


relationship,  Fig.  19,  was 
determined  for  the  lubricating 
oil  from  a  laboratory  sample. 

For  a  test  point,  the 
effective  viscosity  was  found 
from  the  average  temperature 
on  the  bearing  shoe  surface. 
Location  of  the  thermocouples 
is  shown  in  the  sketch. 
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3.  Viscosity  -  Cont'd. 


Turbine  End: 
Free  End: 


Tavg  '  d5  +  T5  +  T?  +  Tg)  ,4 
Tavg  *  +  X2  +  T3  +  I4)  ,4 


For  the  rigid-supports,  each  represented  an  average 
of  6  values  of  temperature  at  the  given  speed  (runs  21-50) . 

For  the  hydrualic- supports,  each  T.  was  the  average  of 
2  values  (runs  1-10) . 

4.  Transmitted  Force 

In  finding  the  trans¬ 
mitted  force  for  a  given 
unbalance  in  the  rotor,  the 
residual  was  mathematically 
subtracted . 


In  the  sketch,  the 
vector  c  is  the  residual. 

The  introduction  of  an 
unbalance  weight  produces  a 
force  a.  With  the  weight 
placed  180°  from  its  original 
location,  b  is  transmitted. 

The  formula  for  finding  d, 
the  component  due  only  to  the 
unbalance,  is  derived  below  by  assuming  a  linear  system. 

b2  =  c2  +  d2  -  2  cd  cos  B 

a2  “  c2  +  d2  -  2  cd  cos  A 

A  +  B  =  180°;  cos  A  ■  -  cos  B 

b2  ■  c2  +  d2  -  2  cd  cos  B 

a2  g  c2  +  d2  +  2  cd  cos  B 

a2  +  b2  -  2  (c2  +  d2)  ~ 

d  -  (a2  +  b2  -  c2  \  1/2 


2 
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Transmitt  ed  Force  -  Cont'd. 
example:  rigid-supports 


a  b  c 


run  21 

154 

run  41 

169 

run  31 

84 

run  26 

159 

run  46 

179 

run  36 

84 

average 

156 

174 

84 

d  -  (  1562  +  1742 
2 


842  ) 


1/2 


142  lb. 


For  a  1.25  lb. -in.  unbalance,  this  is  the  double  amplitude  of 
the  force  transmitted  in  the  X-direction  at  1500  rpm.  However, 
Fig.  15  displays  single  amplitude  values  for  a  unit  unbalance. 


transmitted  force  “  142 _  ■  56.5  lb. 

(2)(1.25)  in. -lb. 


In  addition,  the  force  in  the  Y-direction  is  averaged  with  that 
in  the  X-direction  to  give  a  data  point  in  Fig.  15. 


example:  hydraulic -supports 

In  the  case  of  the  hydraulic-supports,  the  residual  proved 
small  compared  with  the  experimental  error;  therefore,  it  was 
not  necessary  to  use  the  above  procedure  for  subtracting  the 
residual . 

run  1  189 
run  6  220 

average  204  lb.  (double  amplitude  in  X-direction) 

per  3.11  in-lb. 


transmitted  force 


■  204 

(2) (3.11) 


32.8  lb. 

in- lb. 


This  force  was  averaged  with  that  in  the  Y-direction. 
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Transraitted  Force  -  Cont'd. 


wroo 


3  , 

zi.  Eta 


74.  GES 


F 


For  comparison,  it  was 
desirable  to  calculate  the  un¬ 
attenuated  force,  defined  as 
the  bearing  reaction  due  to  a 
static  force  equal  to  the 
unbalance . 

F  -  (27.218) 

g  (74.625) 

For  wr  “  1.0  lb-in.,  the 
unattenuated  force  is  plotted 
in  Fig.  15  as  a  function  of 
speed. 


5.  Vibration  Displacement 

The  calculations  of  journal  vibration  duplicated  those 
for  the  transmitted  force. 


6 . 


Spring  Value  of  Hydraulic -Supports 

ACTUATOR  ACCUMULATOR 


The  spring  value  of 
the  hydraulic -supports  was 
derived  as  shown  below: 

F  a  pA 


t- 


AV  a  -  Aax 

c 

p 


But,  K 

P 


dF 

dx 


(  dV  dF  v 
dV  dx  dp' 
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6.  Spring  Value  of  Hydraulic-Supports  -  Cont'd. 

K  -  LA — E 
P  V 

The  compressed  air  volume  V  is  determined  from  the 
isothermal  relation, 

v  „  PQ  Vo 
P 

1  A2  2 

K  -  LA— El 
P  poVo 

k  -  1.4  (air) 

2 

A  (actuator  piston  area)  ■  14.2  in 

O 

Vq  (accumulator  capacity)  “  30.8  in 

The  accumulator  charge  pressure  was  measured  before  and 
after  each  series  of  testing  over  the  speed  range,  p  is  the 
average  of  data  taken  during  runs  1-10,  with  the  X  and  Y-axes 
considered  together. 

p  represents  the  average  of  the  hydraulic  line  measure¬ 
ments  made  in  runs  1-10  and  in  both  X  and  Y-directions . 

An  average  spring  value  was  calculated  for  each  bearing: 

Free  End: 

K  -  1L.:A)  •  212 (166 .9) 2  .  5510  lb/in. 

P  (46.2)  (30.8) 

Turbine  End: 


K  =  (1.4)(14.2)2(159.8)2  . 

P  (43.7)  (30.8) 


5360  lb/in. 
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ANALYTICAL  COMPARISON 


The  experimental  results  were  compared  to  analytical 
ones,  which  were  determined  using  the  mathematical  tools  developed 
under  this  contract. 


The  generalized-rotor  program?  outlined  in  Part  I  of  this 
volumes  allowed  the  inclusion  of  the  flexibility  and  damping  of  the 
support  together  with  the  spring  and  damping  properties  of  the 
bearings.  The  analysis  in  Vol.  1  of  this  report  enabled  the 
determination  of  the  bearing  properties  from  the  calculated  values 
of  the  Sommerfeld  numbers  given  in  Fig.  13.  The  spring  and  damping 
values?  presented  in  Table  V,  are  equal  in  the  X  and  Y-directions 
because  of  the  symmetry  of  the  bearing  geometry. 


The  rotor  section  data  are  included  in  Table  IV.  Other 
information  used  in  the  generalized-rotor  program  is  given  below: 


Bearing  housing  weight  (free  end) 
Bearing  housing  weight  (turbine  end) 
Support  spring  value  (free  end) 
Support  spring  value  (turbine  end) 
Modulus  of  elasticity 


-  493  lb. 

-  517  lb. 

-  5510  lb/in. 

-  5360  lb/in. 

*  30  x  10°  psi 


For  the  rigid-supports,  the  calculated  journal  response 
is  displayed  in  Fig.  20  together  with  the  experimental  points  for 
comparison.  The  theoretical  curve  shows  no  peaking  within  the 
speed  range.  This  contrasts  with  the  experimental  data  in  which 
peaks  appear  at  2200  rpm  and  4200  rpm.  Calculations  showed  that 
these  resonant  peaks  were  due  to  the  flexibility  of  the  support 
structure  of  between  105  and  5  x  105  lb/in„9  although  a  flexibility 
of  10/  lb/in  was  calculated  from  the  physical  dimensions  of  the 
support  plates.  When  the  transmitted  force  is  compared  in  Fig.  21, 
better  agreement  is  seen  between  the  calculated  and  experimental 
values . 


For  hydraulic- supports?  the  analytical  correlation  of 
the  journal  response  data  is  presented  in  Fig.  22.  The  test  points 
indicate  a  resonance  condition  at  a  higher  speed  than  predicted. 
This  is  the  result  of  the  friction  remaining  in  the  hydraulic- 
supports.  However,  despite  the  friction,  agreement  is  fairly  good. 
At  the  higher  speeds,  both  curves  level  out  to  a  constant  vibration 
amplitude  indicating  the  tendency  of  a  rotor  to  spin  about  an  axis 
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having  the  least  moment  of  inertia. 
The  calculated  curve  shows 
resonant  peaks  at  260  and  440  rpm 
which  represent  the  translational 
and  conical  modes  respectively 
(see  sketch) . 


Fig.  23  compares  the  trans¬ 
mitted  force  levels  for  the  hydraulic 
supports.  The  difference  is  again 
attributed  to  friction  in  the 
supports.  Calculations  assuming  a 
damping  value  of  50  lb-sec/in.  are 
included  to  show  the  effect  of  even 
a  modest  amount  of  damping  on  the 
transmitted  force.  That  is,  the 
resonant  peaks  occur  at  higher 
speeds,  and  the  force  level 
increases  beyond  the  criticals. 

It  is  also  interesting  to  observe 
that  this  damping  drastically 
reduces  the  resonance  peaks. 


T#ANSLAT/OWAL  MOO£~ 


CO/Z/CAL  mode 
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DISCUSSION 


The  test  verified  the  analytical  results  that  in  general 
flexible- supports  could  produce  a  sizeable  reduction  in  transmitted 
force.  Moreover,  it  substantiated  that  hydraulic -supports  were  a 
feasible  approach  for  obtaining  this  flexibility. 

Successful  operation  of  the  test  rotor  up  to  5000  rpm 
yielded  no  unusual  experiences  or  overwhelming  difficulties.  For 
example,  there  were  no  problems  with  instabilities  or  those  intro¬ 
duced  by  resonances  of  subcomponents .  Although  provisions  were  made 
to  lock-out  the  hydraulic-supports  and  to  introduce  additional  damp¬ 
ing  so  that  the  rotor  could  pass  through  the  critical  speeds,  these 
precautions  proved  unnecessary. 

The  main  failing  of  the  investigation  was  the  inability  to 
accurately  predict  the  force  reduction  obtainable  with  hydraulic- 
supports  since  friction  appeared  to  be  inherent  in  the  supports . 
However,  the  force  reduction  advantage  of  flexible- supports  makes 
it  appealing  for  submarine  application.  They  could  be  applied  to 
propulsion  turbines,  turbine - gener at or  sets,  and  other  rotating 
equipment.  For  a  ship's  prototype  design,  an  overhead  coil-spring 

arrangement  (see  sketch  and 
Ref.  4)  is  recommended,  since 
emphasis  should  be  placed  on 
limiting  friction  for  further 
force  reduction  and  more 
accurate  predictability. 
Mechanical  spring  supports 
could  be  designed  with  little 
rubbing  between  mating  parts . 

In  order  to  save  space, 
it  is  further  advised  that  a 
partial  bearing  be  substitu¬ 
ted  for  the  pivoted- shoes 
since  no  oil-whip  instability 
problem  arose  with  its  appli¬ 
cation  in  Ref.  4.  If 
instability  occurred  in  an 
isolated  case,  damping  could  easily  be  added  to  the  supports. 

In  most  applications,  a  servo-return  system,  such  as 
developed  in  Ref.  4,  would  be  incorporated  into  the  supports  to 
correct  for  force  changes  that  would  produce  excessive  deflections 
of  the  rotor .  This  system  could  be  designed  with  sufficient  rugged¬ 
ness  to  insure  reliability.  Moreover,  even  complete  failure  would 
not  interfere  with  the  normal  machinery  operation. 
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RECOMMENDATIONS 

1.  Apply  a  prototype  of  the  flexible-supports  to  a  high-speed  rotor 
onboard  a  submarine. 

2.  Substitute  mechanical  springs  for  the  hydraulic-supports. 
Mechanical  springs  are  a  simplier  scheme  that  would  reduce 
friction  sources. 

3.  Include  an  automatic  positioning  system  if  the  deflections 
on  the  soft  supports  prove  excessive. 

4.  Substitute  a  partial  bearing  for  the  pivoted-shoe  bearing 
in  order  to  save  space  in  an  area  that  will  be  crammed  with 
the  flexible- support  hardware. 

5.  Analyze  the  shock  problem  in  designing  the  ship's  prototype 
for  possible  inclusion  of  snubbers. 

6.  Redesign  the  hydraulic -supports  if  included  in  a  future 
application.  In  order  to  reduce  friction,  it  is  advised 
that  an  overhead  design  with  flat  diaphragm  actuators  be 
considered. 

7.  Defer  the  testing  of  (1)  a  cylindrical  bearing  and/or 
(2)  elastomer-supports.  In  making  the  proposal  for  the 
present  contract,  it  had  been  advised  that  the  desirability 
of  investigating  the  above  items  would  be  determined  at 
the  conclusion  of  the  present  tests.  Now,  little  advantage 
is  seen  in  this  extension  before  the  application  of  the 
present  knowledge  to  a  prototype  design. 


-24- 


CONCLUSIONS 


1.  Hydraulic-supports  proved  to  be  a  feasible  way  of  isolating 
a  rotor,  reducing  the  force  transmitted  to  the  bedplate 
by  as  much  as  a  factor  of  15  when  compared  with  rigid- 
supports  . 


2.  No  instability  condition  arose  when  the  rotor  operated  to 
5000  rpm. 


3.  There  was  no  difficulty  in  passing  through  the  criticals 
that  were  lowered  by  the  soft  supports . 

4.  Friction  in  the  hydraulic-supports  prevented  accurate 
agreement  with  calculated  levels  of  the  transmitted 
force. 
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TABLE  IV 


ROTOR  SECTION  PROPERTIES 


Station 

Weight 

Length 

Cross  Section 

No. 

Lb. 

In. 

In4 

1 

30.0 

5.31 

138.8 

2 

125.0 

6.62 

554.0 

3 

124.0 

2.25 

1102.0 

4 

382.0 

5.91 

3115.0 

5 

876.0 

7.128 

7580.0 

6 

0 

2.652 

7580.0 

7 

777.0 

7.19 

3350.0 

8 

570.0 

8.57 

1554.0 

9 

453.0 

7.19 

1508.0 

10 

235.0 

5.00 

1262.0 

11 

300.0 

6.62 

2130.0 

12 

313.0 

4.87 

543.0 

13 

82.0 

2.13 

188.0 

14 

32.0 

3.19 

117.8 

15 

36.72 

7.42 

20.0 

16 

28.21 

8.65 

5.476 

17 

90.43 

12.375 

.312 

18 

4.297 

12.375 

.12 

19 

4.297 

12.375 

.312 

20 

78.52 

4.75 

2.0 

21 

10.20 

0 

1.0 

Bearing  Stations 

Unbalance 

Station 

1  15  21  6 


-33- 


TABLE  V 

DAMPING  AND  SPRING  VALUES  OF  BEARING  OIL-FILM 


Speed 

rpm 

Free 

End  _  __ 

Turbine 

End 

5- 

ioic 

lb  /in 

3- 
10  C 

lb- sec /in 

105K 
lb/ in 

3— 

10  C 

lb-sec/in 

200 

20.5 

49.0 

18.4 

43.9 

1000 

10.3 

9.8 

9.2 

8.78 

Rigid 

2000 

9.58 

4.9 

8.43 

4.39 

Support 

3000 

8.90 

3.27 

7.97 

2.93 

4000 

7.87 

2.53 

7.66 

2.19 

5000 

7.87 

2.03 

7.20 

1.81 

6000 

7.53 

1.69 

7  05 

1.51 

200 

17.8 

49.0 

15.9 

43.9 

1000 

9.92 

9.8 

8.89 

8.78 

2000 

8.90 

4.9 

8.28 

4.39 

Hydraulic 

3000 

7.87 

3.38 

7.20 

3.02 

Support 

4000 

7.53 

2.53 

6.90 

2.27 

5000 

7.53 

2.03 

6.74 

1.81 

6000 

7.53 

1.60 

6.74 

1.5 

34 


TEST  ARRANGEMENT 


FIGURE  3  BEARING  AND  SUPPORTS 


TEST  ROTOR 
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FIGURE  6  BEARING  PADS  i:  P-69297 
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FIGURE  7 
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HYDRAULIC  SUPPORT 

(looking  fron\  FREE  END') 


Figur-e  <5 


FIGURE  10  RIGID  SUPPORTS 
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INSTRUMENTATION 


FORM  33001  A 


SIGNATURE 


R.J.  Thoman 


FORM  33001  A 


SIGNATURE 


R .  Thoman 


CURVE  NO 


CURVE  NO. 
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Weight  Applied  -  Xb. 


Bearing  Housing  Displacement  -  mils  (.001  In 


SIGNATURE  R‘J«  VW***L 


COMPARISON  OF  FORCE  WITH 


CURVE  NO.. 


H  W 

B  H 
O  t*5 


O 

O 

O'. .  t 

sO 


T> 

«) 

<u 

O  CU 
©  01 
O 

cn  u 
o 
u 

& 


V  r  ^  aaiUITdM^  3T9NIS 

i  -  !:  i  _ .  ;  . . 


FORM  33001  A 


SIGNATURE 


R .  Thoman 


CURVE  NO. 
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NOMRNCIATURE 


a,  b,  c,  d 


g 

k 

P 


r 


t 


w 


X 

A 


A,  B 

C 

C 

s 

D 

F 


K 

K 

P 

L 

N 

P 

R 

S 

T 

V 

X,  Y 
oO 


5 


force  components 
gravitational  constant 
adiabatic  gas  constant 
pressure 

balance  plane  radius 
time 

balance  plug  weight 
displacement 
actuator  piston  area 
angles 

radial  bearing  clearance 
bearing  damping  coefficient 

pedestal  or  support  damping 

bearing  diameter 
transmitted  force 

piston  force 

bearing  spring  coefficient 
pedestal  or  support  spring  constant 

bearing  axial  length 
speed 

bearing  mean  projected  area  pressure 
bearing  reaction 

(F>  **  ;  Sommerfeld  No.,  cycles 

temperature 

volume 

coordinates 

angular  velocity 

viscosity 

phase  angle 


